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ABSTRACT 

We present the first detection of the mean flux of the optical extragalactic background light (EBL) at 
3000, 5500, and 8000A. Diffuse foreground flux at these wavelengths comes from terrestrial airglow, dust- 
scattered sunlight (zodiacal light), and dust-scattered Galactic starlight (diffuse Galactic light). We have 
avoided the brightest of these, terrestrial airglow, by measuring the absolute surface brightness of the 
night sky from above the Earth’s atmosphere using the Wide Field Planetary Camera2 (WFPC2) and 
Faint Object Spectrograph (FOS), both aboard the Hubble Space Telescope (HST). On the ground, we 
have used the duPont 2.5 m Telescope at Las Campanas Observatory (LCO) to obtain contemporaneous 
spectrophotometry of “blank” sky in the HST field of view to measure and then subtract foreground 
zodiacal light from the HST observations. We have minimized the diffuse Galactic light in advance 
by selecting the HST target field along a line of sight with low Galactic dust column density, and 
then estimated the low-level Galactic foreground using a simple scattering model and the observed 
correlation between thermal, 100/im emission and optical scattered flux from the same dust. In this 
paper, we describe the coordinated LCO/HST program and the HST observations and data reduction, 
and present the resulting measurements of the EBL. 

Galaxies brighter than V = 23 AB mag are not well sampled in an image the size of the WFPC2 
field of view. We have therefore measured the EBL from unresolved and resolved galaxies fainter than 
V = 23 AB mag by masking out brighter galaxies. We write as EBL23 to emphasize this bright magnitude 
cut-off. From absolute surface photometry using WFPC2 and ground based spectroscopy, we find mean 
values for the EBL23 of 4.0 (±2.5), 2.7 (±1.4), and 2.2 (±1.0) in units of x 10 _9 ergs s _1 cm -2 sr _1 A -1 in 
the F300W, F555W, and F814W bandpasses, respectively. The errors quoted are la combined statistical 
and systematic uncertainties. The total flux measured in resolved galaxies with V > 23ABmag by 
standard photometric methods is roughly 15% of the EBL23 flux in each band. We have also developed 
a new method of source photometry, uniquely suited to these data, with which we can measure the 
ensemble flux from detectable sources much more accurately than is possible with standard methods for 
faint galaxy photometry. Using this method, we have quantified systematic biases affecting standard 
galaxy photometry, which prevent light from being recovered in isophotes within a few percent of the 
sky level. These biases have a significant effect on faint galaxy counts. The flux from resolved sources as 
measured by our ensemble photometry method is 3.2 (± 0.22), 0.89 (± 0.01), and 0.76 (± 0.01) in units 
of xlO -9 ergs s _1 cm -2 sr _1 A -1 in the F300W, F555W, and F814W bandpasses, respectively, with la- 
combined errors. These values, the total flux from resolved sources, represent absolute minima for the 
EBL23 in each band, and are roughly 30% of the mean flux we measure for the total EBL23. 

Subject headings: Diffuse radiation — cosmology: observations — galaxies: photometry — techniques: 
photometric — interplanetary medium — ISM: dust 


1. INTRODUCTION 

The Extragalactic Background Light (EBL) is the in¬ 
tegrated, mean surface brightness of both resolved and 
unresolved extragalactic sources. In wavelength range of 
these observations (2500-9500A), the EBL is dominated 
by flux from stellar nucleosynthesis at redshifts z < 9, 
and also includes flux from gravitational potential energy 
such as accreting black holes (Fabian & Iwasawa 1999), 
brown dwarfs, gravitationally collapsing systems, and pos¬ 
sibly decaying particles. As such, the EBL is the fossil 
record of the star formation history of the universe and 
a fundamental measure of the luminous energy content of 


the universe. As we show in Figure [|, upper limits from 
previous attempts to measure the EBL and lower limits 
from integrated galaxy counts constrain the EBL to an 
expected level of roughly 1 x 10 -9 ergs s _1 cm -2 sr _1 
A" 1 at 5500A, or 28.2 Ab mag arcsec 2 .[j Because the 
combined flux from foreground airglow, zodiacal light, and 
diffuse galactic light (also plotted in Figure |i) is at least 
100 times brighter than the EBL, a detection of the EBL 
requires measurement accuracies of better than 1%. 

1 A11 surface brightnesses are specified in ergs s —1 cm -2 sr —1 A -1 
unless specifically noted otherwise. AB mag is defined as AB mag= 
—2.5 log F u — 48.6, as usual, with given in ergs s —1 cm -2 Hz -1 . 
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Fig. 1.— Relative surface brightnesses of foreground sources, up¬ 
per limits on the EBL23 (see §[y), and lower limits based on the 
integrated flux from resolved galaxies (V 555 > 23 AB mag) in the 
HDF (Williams etal. 1996). The spectral shape and mean flux of 
zodiacal light and of diffuse galactic light (DGL) are shown at the 
levels we detect in this work. The airglow spectrum is taken from 
Broadfoot &; Kendall (1968) and is scaled to the flux level we ob¬ 
serve at 3800-5100A (see The effective bandpasses for our HST 
observations are indicated at the bottom of the plot. 


Previous attempts to measure the optical EBL have em- 
ployed a variety of different approaches. Mattila (1976) 
attempted to isolate the EBL by differencing integrations 
on and off the line of sight to “dark” clouds at high Galac¬ 
tic latitude, under the assumption that the clouds acted as 
a “blank screens,” spatially isolating all foreground contri¬ 
butions from the background. This pioneering work pro¬ 
duced upper limits and identified both the rapid tempo¬ 
ral variability of terrestrial airglow and extinction and the 
spatial variability of diffuse Galactic light, which proved 
to be the primary obstacles to these early efforts to mea¬ 
sure the EBL. Toller (1983) later attempted to avoid both 
atmospheric and zodiacal foregrounds by using data taken 
with the Pioneer 10 spacecraft at a distance of 3 AU from 
the Sun, beyond the zodiacal dust cloud. Poor spatial res¬ 
olution (2°), however, prevented the accurate subtraction 
of discrete Galactic stars, let alone the diffuse Galactic 
component from these data. Dube, Wilkes, & Wilkinson 
(1977, 1979) made the first effort to measure and subtract 
foreground contributions explicitly based on geometrical 
modeling of airglow and Galactic foregrounds but includ¬ 
ing spectroscopic measurement of the zodiacal light (ZL) 
flux by a technique similar to that which we have adopted. 
Rapid variability caused uncertainty in their airglow sub¬ 
traction which dominated the errors in their results. 

In this work, we take advantage of the significant gains 
in technology and in understanding of the foreground 
sources which have been achieved since the last attempts 
to measure the optical EBL (see Mattila 1990 for a review). 


The most significant technological advance is panoramic, 
linear CCD detectors. Those on-board HST allowed us 
to completely avoid bright, time-variable airglow and pro¬ 
vide sub-arcsecond spatial resolution. High spatial reso¬ 
lution allowed us to resolve stars to V ~ 27.5 mag and 
thereby eliminate the possibility of significant contamina¬ 
tion from unidentified Galactic stars in the held. Ground- 
based spectrophotometry with CCDs also made possible 
much more accurate measurement of the foreground zo¬ 
diacal light than could be achieved with narrow-band Li¬ 
ters and photometers, as were used by Dube, Wickes, & 
Wilkinson (1977, 1979). Finally, IRAS has provided maps 
of the thermal emission from dust at high Galactic lati¬ 
tudes. We have use the IRAS maps to select a line of sight 
for these observations which has a low column density of 
Galactic dust in order to minimize the DGL contribution 
caused by dust-scattered starlight, and also to estimate 
the low-level DGL which cannot be avoided. 

Our measurement of the EBL utilizes three indepen¬ 
dent data sets. Two of these are from HST: (1) images 
from the Wide Field Planetary Camera 2 (WFPC2) us¬ 
ing the F300W, F555W, and F814W Liters, each roughly 
1000A wide with central wavelengths of 3000, 5500, and 
8000A, respectively; and (2) low-resolution spectra (300A 
per resolution element) from the Faint Object Spectro¬ 
graph (FOS) covering 3900-7000A. The FOS data were 
taken in parallel observing mode with the WFPC2 ob¬ 
servations. While Lux calibration of WFPC2 images and 
FOS spectra achieve roughly the same accuracy for point 
source observations, the increase in spatial resolution, 10 4 
times larger held of view, lower instrumental background, 
and absolute surface brightness calibration achievable with 
WFPC2 make it better suited than FOS to an absolute 
surface brightness measurement of the EBL. Nonetheless, 
the FOS observations do provide a second, independent 
measurement of the total background Lux of the night 
sky, also free of terrestrial airglow and extinction, but with 
greater spectral resolution than the WFPC2 images. The 
third data set consists of long-slit spectrophotometry of 
a region of “blank” sky within the WFPC2 held of view. 
These data were obtained at the 2.5m duPont telescope 
at Las Campanas Observatory (LCO) using the Boiler & 
Chivens spectrograph simultaneously with one visit of the 
HST observations (6 of 18 orbits). 

The Low chart in Figure ^ summarizes the reduction 
procedures for each data set used in this measurement, 
the results obtained from each data set individually, and 
the coordination of those results to produce a measure¬ 
ment of the EBL. In this paper, we begin by describing the 
foreground sources in and the details of HST schedul¬ 
ing in (fl The observations and data reduction of both 
HST data sets, WFPC2 and FOS, are discussed in detail 
in this paper. WFPC2 observations and data reduction 
are discussed in §[|. In §[| we present the Lrst results from 
the WFPC2 data, which are measurements of the total 
sky Lux (foregrounds plus background) in each bandpass. 
The FOS observations, data reduction, and results are dis¬ 
cussed in and §[i]. The modeling of diLuse Galactic light 
is discussed in §||. The LCO data and measurement of ZL 
are discussed in Bernstein, Freedman & Madore (2001b, 
henceforth Paper II). A summary of that work is given in 
In §[w[ we present a measurement of the minimum 









Bernstein, Freedman, & Madore 


3 


flux of the EBL from resolved sources, which can be made 
using the WFPC2 data alone. The implications of that 
result are also discussed in Finally, in §|Tl] we com¬ 
bine the results of the individual data sets and modeling 
(horizontal connections shown in the flow chart) to ob¬ 
tain a measurement of the EBL. The implications of these 
results are discussed in Bernstein, Freedman & Madore 
(2001c, henceforth Paper III). 

2. FOREGROUNDS 

The contribution from foreground sources to the flux of 
the night sky is at least two orders of magnitude brighter 
than the expected extragalactic signal. That is, the com¬ 
bined flux from foreground sources is roughly 100 xlO" 9 
ergs s" 1 cm -2 sr" 1 A -1 (or 23.2 AB magarcsec" 2 ) (see 
Figure |l|). A detection of the EBL therefore requires the 
measurement of both the total background and the indi¬ 
vidual foreground sources with an accuracy of roughly 1% 
of the total sky surface brightness, or ~ lxlO" 9 ergs s" 1 
cm -2 sr" 1 A" 1 at 5500A (28.2 AB mag arcsec" 2 ). Each 
of the foreground sources, and our approach to minimizing 
and measuring them, are described briefly in the following 
sections. 


2.1. Terrestrial Airglow and Extinction 

Atmospheric emission is the brightest component of the 
night sky as seen from the Earth’s surface. The mean flux 
and line strengths of the molecular, atomic, and contin¬ 
uum emission which produce “airglow” can vary by several 
percent on time-scales of minutes due to changes in atmo¬ 
spheric column densities throughout the night, stimulation 
by small meteorites, and photo-chemical excitation and 
de-excitation within several hours of twilight (see Cham¬ 
berlain 1966 for background, Dube et al. 1979, and Paper 
II). The mean airglow flux can be predicted based on the 
line-of-sight path-length through the atmosphere, but the 
accuracy of this method is limited to a few percent. Air¬ 
glow subtraction has dominated the errors in several pre¬ 
vious attempts to measure the EBL (see Dube et al. 1979 
and the discussion in Mattila 1990). However, the only 
significant airglow emission seen in the upper atmosphere 
from HST occurs at wavelengths shorter than 2000A and 
only on the daytime side of the orbit (see Lyons et al. 1993 
and reference therein). We have therefore entirely avoided 
optical airglow by using HST as the primary instrument 
in measuring the EBL. 

2.2. Zodiacal Light (ZL) 

Zodiacal light is sunlight scattered off of dust grains in 
the solar system, and it can be as bright as 1500x 10" 9 ergs 
s" 1 cm" 2 sr" 1 A" 1 in the ecliptic plane. The ZL is faintest 
at viewing angles (heliocentric longitude) 130-170 degrees 
away from the Sun, where scattering angles are large, and 
at ecliptic latitudes greater than 30 degrees, where the in¬ 
terplanetary dust (IPD) column densities are lowest. To 
within 10%, the mean flux of the ZL is predictable as a 
function of scattering geometry and ecliptic latitude; how¬ 
ever, this is not accurate enough for our purposes. Instead, 
we have measured the ZL flux directly using detailed a 
priori knowledge of its intrinsic spectrum. Observations 
of the ZL (Leinert etal. 1981, Murthy etal. 1990, Mat- 
suura et al. 1995) show that it has a nearly solar spectrum 


from the UV to the near-IR (1500A-10//m): scattering 
strength is only weakly dependent on wavelength, so that 
the absorption lines and narrow-band features of the Sun 
are very accurately reproduced in the ZL. The broad-band 
color of the ZL is usually described as the ratio of the zo¬ 
diacal to the solar spectrum as a function of wavelength, 


C(X, A 0 ) 


I ZL (X)/I Q (X) 
Izl(Xo)/I®{X 0 ) ’ 


( 1 ) 


in which the reference wavelength, Ao, is typically around 
5500A. Observations of the ZL to date find C(A, Ao) chang¬ 
ing by only 5% per 1000A at the ecliptic orientation and 
scattering angles of interest to us. 

The solar-type spectrum of the ZL is readily explained 
by Mie scattering theory, which describes the scattering of 
light off solid particles larger than the wavelength of the 
incident light. It predicts the mild wavelength dependence 
for scattering by particles with the composition and size 
distribution of the interplanetary dust (IPD) cloud (see 
Roser & Staude 1978, Berriman etal. 1994, and Leinert 
etal. 1998 for a recent summary). Because the equivalent 
widths of the solar Fraunhofer lines are known with high 
accuracy, we can uniquely determine the continuum level 
(mean flux) of the ZL at a given wavelength using the 
apparent equivalent width of the Fraunhofer lines in the 
ZL spectrum. However, high signal-to-noise spectra at 
~ lA resolution are required to measure Fraunhofer lines 
in the spectrum of the ZL. At the time of these observa¬ 
tions, such data could only be obtained from the ground. 
We therefore measured the absolute flux of the ZL using 
spectrophotometry between 3900-5100A obtained at Las 
Campanas observatory simultaneously with the November 
1995 HST observations. To identify the ZL flux over the 
full wavelength range of the HST observations, we com¬ 
bine these observations with an estimate of the ZL color 
obtained from the HST and LCO data together, as de¬ 
scribed in §^| The measurement of the ZL is presented in 
full in Paper II. In §[y], we discuss the use of those results 
in the EBL detection we present here. 


2.3. Discrete Stars and Diffuse Galactic Light (DGL) 

The Galactic contribution to the diffuse night sky comes 
from discrete stars, starlight scattered off interstellar dust, 
and line emission from the warm interstellar medium. 
Light from resolved stars near the optical axis of the tele¬ 
scope can also scatter into the field of view. 

In some previous attempts to measure the optical EBL, 
even light from discrete, resolved stars was difficult to sub¬ 
tract due to poor detector resolution (e.g. Dube etal. 
1979, Toller etal. 1983). In HST/WFPC2 images, how¬ 
ever, Galactic stars in the field can be easily identified and 
subtracted, with negligible residual contribution to the er¬ 
rors in our final results. 

Guhathakurta & Tyson (1989, hereafter GT89), Murthy 
etal. (1990), and others have demonstrated that thermal 
IR emission from Galactic dust and the scattered optical 
flux are well correlated, as both are a function of the col¬ 
umn density of dust and the strength of the interstellar 
radiation field which both heats the dust and is scattered 
by it. We therefore used the IRAS 100/iin maps to select 
a field with low 100/irn emission to ensure that the diffuse 
galactic light (DGL) contribution would be minimal. The 
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Fig. 2.— Flow chart of data reduction, analysis, and results from all data sets used to measure the EBL. Where appropriate, subscripts 
to /(A) indicate the bright magnitude cut-off applied. The thick, horizontal bars divide the pre-reduction and analysis steps. The symbols * 
and ® marking steps in data reduction indicate that a systematic or statistical uncertainty is accrued at that step. Normal type—face indicates 
STScI pipeline data reduction procedures; bold type-face indicates original procedures which we developed for this work. Dashed boxes mark 
estimates of the ZL color; shaded boxes indicate results derived from one data set alone; thick-lined boxes indicate final measurements of the 
EBL from combined WFPC2/LCO and FOS/LCO data sets. 
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Fig. 3.— IRAS 100/im map (Beichmann etal. 1986) covering 
1.5° xl.5°centered on the WFPC2 EBL field. The grey-scale is 
linear: white indicates 0.6 MJy sr -1 at 100 pm, corresponding to 
N(Hl)~ 0.6xl0 20 and E(B — F)~ 0.01 mag; black indicate 1.6 MJy 
sr -1 , corresponding to N(Hl)~ 1.6xl0 20 , E(B — V)~ 0.03 mag (see 
Berriman etal. 1994). HST/WFPC2 and FOS footprints are over¬ 
laid. Stars with 7 < V < 12 mag are marked with circles whose radii 
are linearly proportional the magnitude of the stars. 


approximate 100/im flux in the surroundingl.5° x 1.5° re¬ 
gion of our observations is shown in Figure pi in which the 
position of our target field is shown by the HST WFPC2 
and FOS footprints overlayed the IRAS 100/im map. We 
have estimated the low-level DGL (which contributes even 
along the most favorable lines of sight) by using a simple 
scattering model which accurately reproduces the observed 
scaling relations between 100/im thermal emission and the 
optical scattered light (see Equation [|) . The model com¬ 
bines an estimate of the dust column density based on 
100/im emission with empirical values for the interstellar 
radiation field and dust scattering parameters. Unlike the 
ZL, the flux from DGL is only roughly equal to the EBL 
in surface brightness, so that uncertainty in this modeling 
does not prohibit a detection of the EBL. 

3. FIELD SELECTION AND HST SCHEDULING 

We chose the HST/WFPC2 field (see Table [I| for co¬ 
ordinates) at an ecliptic latitude \(3\ > 30° (to minimize 
contributions from ZL) and near the Galactic pole in a 
region of low Galactic 100/tm emission (to minimize dif¬ 
fuse Galactic light). We also selected the field to avoid 
bright stars, in an effort to minimize the scattered light 
from stars near the optical axis of the telescope. The posi¬ 
tions of stars relative to our field can be seen in the r-band 
CCD image shown in HST Figure || and in the IRAS map 
in Figure 0 . No stars brighter than 12 AB mag fall within 
10 arcmin of the center of our field, and no stars brighter 
than 7 AB mag fall within 2.5 degrees. Based on measure¬ 
ments of the point spread function (PSF) of HST done 
with WFPC1, the attenuation factor at 1 arcmin from the 
center of a point source is 10” 8 (STScI Technical Memos 
RSB-85-03, RSB-85-02, ISR/OTA 06.1). Furthermore, 


• ♦' * 


Fig. 4. — The HST/WFPC2 and FOS footprints overlaid on a 
mosaic Gunn-r—band image (0.6° X 0.6°) taken with the lm Swope 
Telescope at the Las Campanas Observatory. 


the “large-angle” (> 3 arcsec) scattering is roughly one 
order of magnitude lower for WFPC2 than for WFPC1. 
Based on these results, the total contribution from stars 
closer than 5 arcmin to the field is at most 10” 9 of their 
total flux, which amounts to one 35 mag star on-axis. The 
EBL signal, by comparison, is approximately 27.5 mag 
arcsec” 2 . The flux from off axis sources is therefore an in¬ 
significant contribution to the background and a negligible 
source of error. 

The exact scheduling of exposures was also critical to 
this program, as stray solar and terrestrial light in HST 
observations are a strong function of the orbital position of 
the satellite. Sunlight scattered off the limb of the Earth 
can increase the background level by a factor of 10 when 
the bright limb is near the optical axis of the telescope. 
All of our science observations were therefore scheduled to 
execute exclusively in the shadow of the Earth. Because 
our field was at a viewing angle angle greater than 135° 
from the Sun during the months of our observations, the 
satellite was pointed in a direction greater than 45° away 
from the satellite’s direction of motion when our field was 
observed from within the Earth’s shadow. This guaran¬ 
teed that no upper-atmosphere glow from the satellite’s 
flight path would affect our observations. Scattered moon¬ 
light was also explicitly avoided by exposing only with the 
Moon greater than 65° from the optical axis of the tele¬ 
scope, which is the angular separation at which models 
predict that the attenuation function for off -axis scattered 
light becomes flat (STScI Technical Memo RSB-85-03). 
Based on those models, confirmed by on-orbit measure¬ 
ments (Burrows 1991, Hasan & Burrows 1993), the light 
from the moon at 65° from the optical telescope assembly 
(OTA) is less than 10” 5 photons s” 1 arcsec 2 per 100A, 
which is 1000 times smaller than the lower limit expected 
for the EBL, and is therefore also a negligible source of 
background during our observations. 
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Table 1 

EBL Field Coordinates 


System 

longitude 

latitude 

Equatorial(J2000) 

3 h 00 m 20?40 

-20°10'45'/3 

Galactic 

206? 6 

-59? 8 

Ecliptic 

35? 5 

-35? 5 


Our 18 orbits were divided into three visits, staggered 
by one month each between October and December, 1995. 
The ZL flux changes with the Earth’s orbital position 
about the Sun and the line of sight through the interplan¬ 
etary dust during the three visits. We therefore expected 
to see changes in the total sky flux between visits, at a 
level predicted by the models of the interplanetary dust. 
During all three months, our target held has a ZL surface 
brightness which is within 20% of the absolute minimum 
intensity of the ZL at any orientation. By separating the 
orbits for this program into three visits, we were also able 
to observe the held at two different roll angles - V3 posi¬ 
tion angle 222° for the November and December visits, V3 
angle 132° for the October visit — thereby changing the 
position of the off-axis stars with respect to the OTA. This 
allows us to further rule out signihcant off-axis scattered 
light and other possible photometric anomalies. 

4. HST/WFPC2 DATA 
4.1. Observations 

The 18 orbits scheduled for this program executed on 
28 October, 29 November, and 16 December, 1995. Dur¬ 
ing each orbit, a single long integration (1800 sec) was 
taken while the satellite was in the Earth’s shadow and 
one additional short integration (300-400 sec) was taken 
to hll the time between target acquisition and entry into 
the Earth’s shadow. Two orbits per visit were spent ob¬ 
serving with each of three WFPC2 filters (F300W, F555W, 
and F814W). The short exposures were combined and used 
to improve cosmic ray rejection. In all, five 1700 sec darks 
and ten bias images were taken during the bright portion 
of the six orbits comprising each visit. The gain setting 
for these data was 7 e _ /DN (bay 4) and the read noise was 
roughly 5.2 e~ for each of the four CCDs. The field of view 
of each WFPC2 image is roughly 4.4arcmin 2 (725 x 725 
well-exposed pixels in each of the 3 WF CCDs, with 0.0996 
arcsec 2 per pixel). 

4.2. Basic Data Reduction and Calibration 

The reduction of these data involved a mixture of the 
standard pipeline and our own methods. The pipeline pro¬ 
cedures used included the corrections for analog-to-digital 
conversion errors, overscan subtraction, bias image sub¬ 
traction, flat-fielding, and point-source calibration. The 
accuracy of these steps is largely provided in STScI doc¬ 
umentation. Whenever possible, we have independently 
verified the errors we quote. All steps in the reduction 
are indicated in the flow chart in Figure with associ¬ 
ated errors explicitly included in our error budget. The 
reduction is discussed briefly below; greater detail can be 
found in Bernstein (1998). The following section should be 


read by those who are interested in the reduction proce¬ 
dures and detailed error discussions. As the goal of these 
observations is an absolute measurement of the total flux 
of the night sky, zero point calibration is as critical to 
the accuracy of the results as flux calibration. Table || 
summarizes the errors discussed below and indicates the 
dominant sources of error. 

4.2.1. Pipeline 

The WFPC2 analog-to-digital conversion (ADC) has a 
bias towards DN values which correspond to the setting 
of all low-order bits in the digitization, with larger sig¬ 
nal levels being more strongly affected. The correction for 
this bias involves replacing the output value with a value 
corresponding to the scaled average input signal for that 
output. DN values in the range of our data have an av¬ 
erage additive correction of 0.86 DN, with the correction 
smoothly varying for levels in this range. The same cor¬ 
rection is used for each of the four CCDs in WFPC2. The 
data used to determine this correction were taken pre- and 
post-dynamic testing, with only 0.02 DN variation. The 
correction is thought to be stable under normal usage and 
insensitive to temperature fluctuations (STScI Technical 
Memo RSB-85-01). At the signal levels of our data, the 
error after the ADC correction applied is 0.02 DN. 

The bias level was removed with a two-step process of 
overscan subtraction and bias-image subtraction. We em¬ 
ployed the usual method established for WFPC2 data of 
subtracting overscan from odd and even columns using the 
overscan from each individual exposure. For bias-image 
subtraction, we used the “superbias” frames produced 
for the reduction of the Hubble Deep Field (HDF, see 
Williams etal. 1996). This bias frame was produced from 
the average of 200 frames, and has lower read noise than 
we could achieve using the 30 bias frames taken during 
our orbits. Instead, bias frames from our orbits were test- 
reduced (overscan and bias subtraction using the HDF su¬ 
perbias) to verify that the superbias produced corrected 
levels consistent with zero, with a scatter of 0.002 DN. As 
this error is inseparable from our estimate of the error in 
the dark subtraction, it is included with the dark subtrac¬ 
tion in our error budget. 

We used the pipeline flat-fielding images provided by 
the WFPC2 GTO team in May 1996. The pixel-to-pixel 
flux error in those images is reported in the WFPC2 Data 
Handbook (V3.0) to be roughly 0.3% ( rms ) for the WF 
chips, and 0.5% for the PCI. The errors over spatial scales 
greater than 10 arcsec are less than 0.5% for all four chips. 
As we excluded data within 75 pixels of the edge of each 
chip from our analysis, issues of geometrical distortion and 
vignetting are avoided. The photometric calibration of 
WFPC2 is tied to the mean level of the flat- fielding im- 
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Table 2 

Sample of V 555 Source Catalog 


ID# 

X 

y 

RA(h:m:s) 

Dec(h:m:s) 

mtot 


^■iso 

^(^iso) 

A . 

-'MSO 

a 

elong. 

EBL555-2-1 

129 

55 

3:00:21.66 

-20:10:37.4 

25.51 

0.04 

25.45 

0.05 

11 

0.53 

1.08 

EBL555-2-2 

483 

55 

3:00:24.14 

-20:10:35.4 

25.94 

0.10 

26.06 

0.09 

20 

1.72 

2.04 

EBL555-2-3 

623 

58 

3:00:25.13 

-20:10:35.1 

26.32 

0.13 

26.37 

0.10 

16 

1.50 

2.03 

EBL555-2-4 

74 

60 

3:00:21.28 

-20:10:38.2 

27.79 

0.26 

28.27 

0.31 

3 

0.68 

1.57 

EBL555-2-5 

600 

61 

3:00:24.97 

-20:10:35.4 

26.62 

0.13 

26.71 

0.12 

11 

1.04 

1.43 

EBL555-2-6 

279 

63 

3:00:22.72 

-20:10:37.3 

27.12 

0.17 

27.09 

0.15 

9 

0.83 

1.15 

EBL555-2-7 

298 

57 

3:00:22.85 

-20:10:36.6 

25.46 

0.20 

26.16 

0.11 

29 

4.81 

7.60 

EBL555-2-8 

306 

65 

3:00:22.91 

-20:10:37.3 

27.47 

0.24 

27.73 

0.23 

6 

1.08 

3.65 

EBL555-2-9 

155 

68 

3:00:21.85 

-20:10:38.5 

26.60 

0.12 

26.70 

0.12 

10 

0.95 

1.29 

EBL555-2-10 

302 

69 

3:00:22.88 

-20:10:37.7 

27.22 

0.34 

27.70 

0.23 

7 

1.15 

2.62 

EBL555-2-11 

650 

68 

3:00:25.32 

-20:10:35.9 

25.09 

0.05 

25.08 

0.05 

33 

1.30 

1.09 

EBL555-2-12 

756 

70 

3:00:26.06 

-20:10:35.7 

25.60 

0.12 

25.93 

0.10 

34 

2.33 

1.53 

EBL555-2-13 

380 

75 

3:00:23.43 

-20:10:37.9 

26.75 

0.12 

26.70 

0.11 

9 

0.79 

1.25 

EBL555-2-14 

393 

74 

3:00:23.53 

-20:10:37.8 

27.42 

0.21 

27.40 

0.17 

7 

0.82 

1.39 

EBL555-2-15 

729 

75 

3:00:25.88 

-20:10:36.3 

27.24 

0.20 

27.51 

0.20 

7 

0.92 

1.67 

EBL555-2-16 

680 

77 

3:00:25.54 

-20:10:36.7 

26.87 

0.17 

27.29 

0.17 

9 

1.15 

1.85 

EBL555-2-17 

759 

82 

3:00:26.09 

-20:10:36.8 

27.79 

0.49 

28.22 

0.33 

4 

0.95 

1.76 

EBL555-2-18 

534 

85 

3:00:24.52 

-20:10:38.1 

25.83 

0.12 

26.22 

0.11 

28 

1.79 

1.31 

EBL555-2-19 

86 

89 

3:00:21.37 

-20:10:40.9 

27.05 

0.19 

27.22 

0.16 

8 

0.93 

1.41 

EBL555-2-20 

178 

95 

3:00:22.02 

-20:10:41.0 

26.72 

0.17 

26.97 

0.15 

10 

1.15 

1.35 

EBL555-2-21 

96 

98 

3:00:21.45 

-20:10:41.8 

27.56 

0.23 

28.02 

0.28 

4 

0.72 

1.45 

EBL555-2-22 

670 

100 

3:00:25.48 

-20:10:39.0 

27.19 

0.23 

27.56 

0.22 

10 

0.94 

1.21 

EBL555-2-23 

321 

100 

3:00:23.03 

-20:10:40.7 

26.30 

0.12 

26.76 

0.14 

13 

1.16 

1.51 

EBL555-2-24 

276 

105 

3:00:22.71 

-20:10:41.4 

27.72 

0.28 

27.81 

0.21 

6 

0.75 

1.71 

EBL555-2-25 

599 

106 

3:00:24.98 

-20:10:39.8 

27.07 

0.21 

27.31 

0.18 

10 

1.02 

1.40 

EBL555-2-26 

213 

107 

3:00:22.27 

-20:10:42.0 

25.91 

0.13 

26.32 

0.11 

24 

2.19 

1.93 

EBL555-2-27 

666 

112 

3:00:25.45 

-20:10:40.1 

26.98 

0.21 

27.31 

0.19 

9 

1.32 

3.09 

EBL555-2-28 

311 

113 

3:00:22.96 

-20:10:42.1 

27.48 

0.23 

27.80 

0.21 

6 

0.94 

2.00 

EBL555-2-29 

645 

106 

3:00:25.30 

-20:10:39.6 

22.01 

0.01 

22.01 

0.01 

375 

4.82 

1.46 


Note: Coordinates are J2000. Total magnitudes (m to t), isophotal magnitudes (TOi so ), and the associated rms 
errors (a) are in AB magnitudes. The isophotal area (Aj so ) and flux-weighted, second order moment along the 
major axis (a) are in units of WF pixels, which are 0.0996 arcsec on a side. The isophotal, weighted elongation 
(e) is the ratio of the second order moments along the major and minor axes, as calculated by SExtractor. 
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ages, so that no systematic error is introduced due to flat- 
fielding. 

4.2.2. Subtraction of Dark Backgrounds 

We have developed a new method for dark subtrac¬ 
tion because the pipeline subtraction ignores the fact that 
the signal which accumulates on the WF chips while the 
shutter is closed varies between 0.5-1.1 DN per 1800 sec 
exposures for the WF3 and WF4 chips (more for WF1, 
less for WF2). This variability is not caused by ther¬ 
mal dark current, but rather by variations in the dark 
“glow” contributed by the MgF 2 field flattening lenses, 
positioned immediately above each of the WFPC2 chips. 
When struck by even low-energy (200-2000 eV) ions and 
electrons (which are abundant in the upper atmosphere 
in the form of and cosmic rays), MgF 2 will produce 
broad-band fluorescence on a time-scale << 1 sec (Qi 
etal. 1991). We have found that the cumulative flux in 
cosmic ray hits recorded in a dark (shutter-closed) expo¬ 
sure correlates strongly with the mean dark level in that 
frame. This correlation is shown Figure [ 5 ] in which we plot 
the cosmic ray flux and dark glow for each WFPC2 CCD 
in each of 15 dark exposures. The scatter around the linear 
fit for each chip is ~ 0.06 DN. This plot also shows clearly 
that the dark background in a 1700 sec exposure can vary 
by as much as ±0.35 DN. As the total background sky sig¬ 
nal through the F300W is only 0.3 DN/1700sec, accurate 
subtraction of the dark glow is essential for this measure¬ 
ment. 

In collaboration with H. Ferguson (STScI), we devel¬ 
oped a method which uses the correlation between the 
MgF 2 glow and the flux in cosmic rays to isolate the ther¬ 
mal dark component from the dark glow. We obtained a 
final dark solution using 80 individual dark images taken 
between October and December 1995, and have tested the 
accuracy of this method by test-reducing the 15 dark expo¬ 
sures taken between the science exposures of our program. 
These 15 darks were not used in determining the dark so¬ 
lution itself. After ADC correction, overscan subtraction, 
bias subtraction, and dark subtraction by the prescription 
described above, the average mean level of the 15 darks 
is consistent with zero, with a scatter of < 0.05 DN per 
1700 sec exposure. The statistical error in the pipeline 
dark-subtraction, by comparison, is roughly 0.15-0.25 DN 
for the WF3. 

4.2.3. Charge Transfer Efficiency (CTE) 

Accurate photometry depends on stable charge trans¬ 
fer efficiency (CTE), which is the efficiency with which 
electrons are transferred between pixels during read out. 
Variable CTE as a function of total charge level, time, 
or position over the chip results in non-linear sensitiv¬ 
ity. Unfortunately, the WFPC2 chips are known to have 
a “a small parallel CTE problem” (Holtzman et al. 1995b, 
henceforth H95b; Whitmore, Heyer, & Casertano 1999): 
the percentage of electrons which are read-out of a pixel 
depends on the row number of the pixel in question, the 
number of electrons in that pixel when read-out begins, 
and the mean charge level in the pixels through which the 
charge travels during read-out. 

Based on in flight point source photometry and on lab¬ 
oratory tests using a CCD from the same silicon wafer as 



Fig. 5. — Correlation between the cumulative flux in cosmic ray 
events and the mean dark level (thermal dark plus MgF 2 glow) for 
each WFPC2 chip in 15 exposures: x’s for PCI fluxes, open boxes 
for WF2, circles for WF3, and filled circles for WF4. Statistical 
errors in both quantities plotted are negligible; scatter in the cor¬ 
relation probably reflects statistical variations in the energy level of 
the cosmic rays hitting the CCD and the MgF 2 plates. 


WFPC2 CCDs, the CTE variability is known to be caused 
by electron traps which were in the silicon itself before the 
pixel mask was etched into the wafer (J. Trauger, private 
communication). Although electrons can be trapped in 
the silicon at the pixel where they are detected, a greater 
surface area of silicon is encountered during read-out, pro¬ 
viding greater opportunity for trapping if the traps are not 
already filled before readout begins. When a bright point 
source is imaged against a faint background, many new 
traps are encountered during read-out. For example, 4% 
of a 10,000 e~ point source are lost when the charge is 
transferred over 800 rows if the background level is ~10 
e _ /pixel. If the traps are already filled by a high back¬ 
ground level over the whole chip, fewer electrons will be 
lost from sources during read-out. Surface photometry is 
less affected by this CTE problem because all pixels are 
filled to the same charge level and, hence, new traps are 
not encountered during read-out. 

The signal level read out for our images is roughly 80e- 
/pixel for the F555W and F814W images and 2e-/pixel for 
the F300W images. We have identified the CTE losses in 
these data by two methods. First, data from WFPC2 and 
laboratory tests conducted by J. Trauger of point sources 
show that ~0.35 e _ /pixel are lost from a point source 
which is imaged against a background level of roughly 
80 e _ /pixel, while ~l.le _ /pixel are lost when the back¬ 
ground level is zero (J. Trauger, private communication). 
The difference between the number of electrons lost at the 
two different background levels indicates that the number 
of single-electron traps available in a uniform background 
of 80e _ /pixel is 0.75 e _ /pixel (~l.l-0.35). Second, to con¬ 
firm that this trapping reflects the CTE losses for uniform 
sources, we have conducted another set of tests with the 
help of J.Trauger to measure the non-linear response of the 
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Fig. 6. — Histogram of cosmic ray events as a function of the 
number of affected pixels (top panel), total flux (middle panel), and 
peak flux (bottom panel). The number of events falls off in each plot 
before an event would become difficult to detect in these data. 


spare WFPC2 CCD to uniform, low-level backgrounds. 
In this test, we simply expose the CCD to a light source 
which is stable to better than 0.1% for varying lengths 
of time and look for non-linearities in the detection rate. 
The results of these uniform-source test are in excellent 
agreement with results from point-source tests. We there¬ 
fore conclude that 0.75 (±0.25) e _ /pixel are lost when a 
mean level of ~ 80e“/pixel are read-out, and 0.1 (±0.05) 
e _ /pixel are lost from a mean level of ~ 2e _ /pixel. The 
maximum uncertainty in these corrections is a negligible 
contribution to our final results (see Table ||). For more 
detail on the result of the laboratory tests, see Bernstein 
(1998). 

4.2.4. Cosmic Ray Detection 

If low-energy cosmic ray events were left unidentified in 
WFPC2 data, we would unintentionally include the signal 
from those events in the inferred background level. Figure 
|(] shows the total number of cosmic ray events as a function 
of total flux, peak flux, and affected pixels per event. In 
each histogram, the number of cosmic ray events falls off 
sharply well before cosmic rays would become difficult to 
detect: both peak and total fluxes for cosmic ray events are 
significantly higher than the mean level of the background 
(< 10 DN), and the total number of pixels affected by 
a single cosmic ray event rarely drops below three. We 
therefore conclude that the cosmic ray removal process 
detects virtually all cosmic rays which hit the chips. 

4.2.5. Flux Calibration 

The final accuracy of a surface brightness measurement 
depends on four independent facets of the calibration: 
point source calibration; aperture correction, compensat¬ 
ing for flux not recovered in the point source calibration; 
the calibration of the fiducial standard star system; and 
the solid angle subtended by each pixel. The solid an¬ 
gle of the pixels is, of course, very well known and is 


a negligible source of error for the WFPC2 data. How¬ 
ever, the other three calibration steps require quantities 
of data and observing time only available to the WFPC2 
Instrument Team. Details of the WFPC2 calibration are 
thoroughly discussed in Holtzman et al. (1995a, henceforth 
H95a), H95b and subsequent STScI Instrument Science 
Reports. Below, we summarize their results. 

Sensitivity to the same fiducial standard is stable to 
a < 1% for the F555W filter, and a < 1.5% for the 
F814W and F300W between February 1995 and March 
1997 (STScI TIR/WFPC2 97-01). Observations of the 
same standard may be variably affected by CTE due to 
variations in the peak flux of the observation which will 
arise from variations in the position of the PSF within 
a single pixel. The real response of WFPC2 is, there¬ 
fore, likely to be more stable than this result. The syn¬ 
thetic photometry package (SYNPHOT), including system 
throughput curves and synthetic zero-points, is based on 
the HST secondary standard system and is estimated by 
H95b to be accurate for point sources to roughly 1% for a 
source with typical stellar colors. Point-source calibration 
should therefore be good to roughly 1%. 

Aperture corrections based on the encircled energy curves 
are given in H95a (see Figures 5a & 5b, and Tables 2a & 
2b in H95a). For the optical filters, the aperture correction 
is nearly 10% in moving from a 0.5arcsec aperture to an 
infinite aperture, with uncertainties of 1% for the redder 
bandpasses, and 1-2% for the F300W. 

Finally, the secondary standard star system which is 
used for all the HST instruments is “conservatively esti¬ 
mated” to be within 1% of the Hayes 1985 optical cal¬ 
ibration of Vega (Bohlin 1995; Colina & Bohlin 1994; 
Hayes 1985), with internal agreement of the same order. 
The agreement between the HST secondary system and 
the Hamuy etal. (1992) secondary system to which our 
ground-based observations are tied is estimated to be 1- 
2%, as discussed in Paper II. 

4.3. Combined Images 

We produced combined images in which individual re¬ 
solved sources could be identified and photometered for 
two purposes: (1) to allow bright galaxies and stars to be 
identified and excluded from the EBL measurements, (2) 
to allow comparison of the galaxy sample and EBL results 
to other data sets. The measurements of the mean sky 
flux discussed in and § [Tc| are obtained from individ¬ 
ual WFPC2 exposures. The deepest combined image for 
each filter is obtained from the four exposures taken during 
the November and December visits, during which the tele¬ 
scope roll angles were the same. The small differences in 
the alignment of the WF CCDs prevent data taken during 
the October visit from being usefully included in combined 
images for the purposes of accurate photometry. The com¬ 
bined images and resolved source catalogs are available on 
the NASA/IPAC Extragalactic Database (NED).^] A sam¬ 
ple of the catalogs is shown in Table || and the combined 
F555W image is shown in Figure [7] The F300W, F555W, 
and F814W catalogs contain 140, 687, and 644 objects, 
respectively. The photometry of these sources is discussed 
in the next subsection. 


2 Web address: ittp://nedwww.ipac.caltech.edu/level5/ 
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Table 3 

WFPC2 Background Flux: Errors Per Image 



F555W, F814W 

F300W 

Statistical (per image) 

Off-axis scattered light (§g) _ _ 

< 0.001 DN 

< 0.001 DN 

Dark and bias subtraction (§§1.2.2, 4.2.1) 

±0.05 DN 

±0.05 DN 

Flat—ficlflinf (44.2.1) 



CTE (§ 4.2.3) 

±0.03 DN 

±0.007 DN 

Cumulative Statistical Error * 1 

0.06 DN 

0.05 DN 

Systematic 

A-to-D conversion b (§ ip. 1|) 

Point source flux c.a.l.f bk.‘2., ) 

±0.02 DN (0.2%) 

±0.02 DN ( 6 %) 

1 % 

1.5% 

Aperture correction (S4.2 5) 

1 % 

1 % 

Solid angle (pixel scale) (§ 1.2.5) 

0 . 1 % 

0 . 1 % 

Combined Systematic Error 1 ' 

~T7Z% 

5.6% 


individual random errors have been added in quadrature to obtain a cumulative, 
one-sigma random error in DN. To obtain fractional errors, compare to the mean 
flux per exposure (~0.3 DN per pixel for the F300W, ~12.5 DN per pixel for F555W 
and F814W). 

b The A -to-D conversion is an additive correction, the associated error is therefore 
the ratio of the A-to-D error and the mean level in the frame. All other systematic 
uncertainties listed are multiplicative. 

c Systematic errors have been combined assuming a flat probability distribution 
for each contributing source of error. The resulting systematic error is roughly 
Gaussian distributed, and the quoted value is the 68 % confidence interval. For a 
detailed discussion see § |Tl| . 



Fig. 7.— The combined F555W images from November and De¬ 
cember visits. The total exposure time is 4 x 1800 sec. 


4.4. Resolved Objects 

We used the image analysis package SExtractor (Bertin 
& Arnouts 1996) for object identification and photometry. 
The requirement for detecting an object in our data is a 
core surface brightness, /i core , of at least 2 <r s k y within the 


core 4 pixels. The isophotal detection threshold, p,i so was 
defined as 1.5<r s ky over 6 pixels, after smoothing with a 
3x3 boxcar filter. This corresponds to isophotal detec¬ 
tion thresholds of 24.7, 25.8, and 25.3 ST mag arcsec -2 , at 
F300W, F555W, and F814W, respectively. The shape of 
the smoothing kernel and the detection threshold have lit¬ 
tle affect on detection of objects with V 555 < 27.5 AB mag, 
and photometry is not strongly affected for objects with 
V 555 < 27 AB mag. Total magnitudes are defined as the 
flux within an area at least three times larger than the 
ellipse defined by the first moment radius (see Kron 1980) 
and isophotal elongation. For well-detected sources, the 
detection and photometry parameters we employ are sim¬ 
ilar to those used by the HDF team in producing the 
HDF catalogs (see the photometry discussion in Williams 
etal. 1996). Isophotal detection limits are roughly 1.5 mag 
brighter than the HDF, in keeping with the difference in 
exposure time, which is a factor of 15 — 21 in the various 
bandpasses. Differential number counts for our held are 
plotted in Figure |[ Our detection limits are illustrated 
by the turn-over in the counts in each band. The y/N— 
error bars show that V ~ 23 AB mag is reasonable bright- 
magnitude cut-off for observations in a held of this size. 
The galaxy counts from the HDF are also plotted in Figure 
|s| to demonstrate that the held we have observed contains 
a typical held-galaxy population in both color and num¬ 
ber density. The HDF is a convenient benchmark for this 
comparison simply because the data are publicly available, 
widely studied, and of a similarly “blank” held. The HDF 
counts show twice as many galaxies at magnitudes brighter 
than 23 AB mag in each of the three band, while the differ- 
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Fig. 8 .— Differential galaxy counts with y/N error bars from WFPC2 EBL images (filled circles), and the HDF images at F300W, F606W, 
and F814W (x’s). HDF counts are taken directly from the published catalogs (Williams etal. 1996). The difference in completeness limits 
between the two data sets reflects the difference in total exposure times, which are roughly x21, xl 6 , and xl7 at U , V, and I, respectively. 
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Fig. 9.— The flux from bright (V 555 < 23 AB mag) sources 
relative to the total flux in each WFPC2, EBL image is plotted 
for four exposures through three filters (twelve images total). The 
total flux is just the sum of all counts in the calibrated images. 
The flux from bright sources is a background-subtracted flux w ithin 
four times the isophotal radius (4r; ao ) for galaxies (see § [l 0. l| ) and 
within a 5 arcsec radius for stars. The F555W, F814W, and F300W 
flux ratios are marked by open circles (overlapping), triangles, and 
crosses, respectively. 


WF3 
WF CCD 


ence in number density fainter than 23ABmag is modest. 
Such differences at bright magnitudes are well within the 
typical field -to-field fluctuations for galaxy counts and il¬ 
lustrate the need for a bright magnitude cut-off. 

Four stars are detected in the WF chips, three in WF2 
with V 555 = 19.0, 20.5 and 20.8 AB mag, respectively, and 
one in WF4 with V 555 =22.0 AB mag. Star-galaxy sep¬ 
aration poses no difficulty; stars brighter than V 555 ~ 
23.0 AB mag are masked out, regardless, and the flux from 
stars beyond that limit is < 10 % of that from detected 
galaxies at the same apparent magnitudes (Infante 1997). 

4.5. On-Axis Scattered Light 


Scattered light from the “bright” (19 < V < 23 AB 
mag) objects imaged in the field is one possible source of 
error in a measurement of the background from fainter ob¬ 
jects. To quantify this effect, we begin by considering how 
much light these sources produce relative to the total flux 
in an image. In Figure we plot the ratio of the flux from 
bright sources to the total flux detected in each of the three 
WF CCDs during each of four exposures. The total flux 
is just the sum of all counts in the calibrated images. The 
flux from bright sources is a background-subtracted flux 
withi n fou r times the isophotal radius (4ri so ) for galaxies 
(see j jlO.l] ) and within a 5 arcsec radius for stars. For each 
bandpass, the ratio of flux in bright sources to the total 
flux varies from chip to chip. In the WF2 or WF3 chips, 
bright objects contribute only ~ 5% of the total integrated 
flux. As several of the brightest galaxies in the field are 
imaged on the WF4 chip, almost ~ 15% of the flux in 
F555W and F814W in that field comes from objects with 
V < 23 ST mag. 

In order exclude these V < 23 AB mag objects from 
the measured background, we masked regions around each 
resolved galaxy which extend to four times the isophotal 
detection radius of the galaxy, and masked regions with 
radii of 5 arcsec for stars. Based on the encircled energy 
curves in H95a and growth curves in our own images (see 
j) |lo| and Appendix B), we estimate that less than 2% of 
the light from these objects is imaged beyond the extent of 
the masked regions. As an upper limit on their scattered 
light, 2 % of the flux from these objects constitutes only 
O.lx 10 -9 ergs s -1 cm -2 sr _1 A -1 , less than 3% of the total 
flux at 5500A in our estimate of the EBL. As discussed 
previously in §[| the scattered light from off -axis sources 
is also negligible. 

At the wavelength of the F300W filter, the flux from 
bright objects on the WF4 chip is 15-20% of the total 
flux, and variation in the percentage flux contributed by 
objects in the F300W images is greater than in the other 
two bands. These characteristics are explained by two 
facts: ( 1 ) errors in dark glow subtraction cause noticeable 
variation in the background level between exposures in this 
band; and ( 2 ) the ZL falls off rapidly below 4000A, so that 
the ZL contribution to this band is a factor of three smaller 
than in the F555W or F814W and variations in source flux 
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Fig. 10. — Combined F555W WF2 image. North is up and East 
is left. The ellipses indicate the region extending to 4x the isophotal 
detection radius, as discussed in §llCJ. Red ellipses mark galaxies with 
V 555 < 23 AB mag; green ellipses! - !;! < V 555 < 25 AB mag; and blue 
ellipses, 25 < V 555 < 28.5 AB mag. All sources in the catalog are 
shown. The completeness limit is V 555 ~ 27.5. 


from chip to chip are fractionally larger in this band than 
in the other two. In contrast, ZL contributes roughly 95% 
of the diffuse background at 5500-9000A (see §| and Paper 
II). Nonetheless, the contribution to the background from 
detected objects in the F300W images, again assuming 2% 
scattered light, is still a minor uncertainty in our results, 
being 0.3x 10 _9 ergs s -1 cm -2 sr -1 A -1 , or 7% of the EBL 
we detect at 3000A. 

5. HST/WFPC2 TOTAL MEASURED BACKGROUND 

After calibration, we measure the total background sig¬ 
nal by simply taking the average of the detected flux per 
pixel, excluding pixels which fall into any of the following 
three categories: those flagged as bad in the data qual¬ 
ity file for any reason (see WFPC2 Data Handbook V3.0), 
those within two pixels of a cosmic ray event, or those 
within the masked region associated with a star or galaxy 
brighter than V 555 = 23 AB mag, as described in the previ¬ 
ous section. We show the masked regions in Figures [L0] ^2|, 
for each of the WF chips. Masks were defined based only 
on the F555W images; the same object masks were ap¬ 
plied to all three bandpasses. The detected EBL at all 
wavelengths is therefore defined by contributions from the 
same sources. 

Each WF image from each 1800 sec exposure produces a 
measurement of the mean background. We have averaged 
the six measurements from each visit to obtain a single 
background measurement from each visit. The results for 
the November and December 1995 visits are compared in 
Figure |l3] to illustrate the change in background due to 
the geometric path length through the IPD, and thus the 
difference in ZL contribution during the two sets of ob¬ 
servations. The same modulation in background flux was 
identified in the FOS data between November and Decem- 


Fig. 11.— Combined F555W WF3 image. North is up and East 
is left. Ellipses are as in Figure 10. 



Fig. 12.— Combined F555W WF4 image. North is up and East 
is left. Ellipses are as in Figure 10. 


ber 1995, as discussed in (171 The predicted change from 


October to November is oforder 10%. However, as there 
are no FOS or ground-based measurements in October, 
the interpretation of the modulation as a change in ZL is 
not conclusive for that epoch. We do not include the PC 
data in the final average because the dark subtraction is 
least accurate for that chi p due to the smaller pixels and 
higher dark current (see § [4.2. 2 ), and because the field of 
view of the PC is a negligible increase over that of the 
other three chips. However, the PC chip gives the same 
result as the WF chips to within 2%, which is consistent 
with the variation in object distribution over the 4 chips 
and dark subtraction errors. 

The error bars plotted in Figure O show the statistical 
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Table 4 

Total Sky Flux (Foregrounds and EBL) Measured from WFPC2 


Filter 

Ao(FWHM) A 

Total Background 

(± rms) a 

[± sys] b 

November 1995 

F300W 

3000(740) 

3.35x10“® 

(± 9%) 

[± 5.6%] 

F555W 

5500 (1230) 

I.O 6 XIO - 7 

(± 0.3%) 

[± 1.4%] 

F814W 

8000(1490) 

7.24xl0 -8 

(± 0 . 2 %) 

[± 1.4%] 

December 1995 

F300W 

3000(740) 

3.12x10-® 

(± 7%) 

[± 5.6%] 

F555W 

5500(1230) 

1.04xl0 -7 

(± 0.4%) 

[± 1.4%] 

F814W 

8000(1490) 

7.11x10-® 

(± 0 . 2 %) 

[± 1.4%] 


a Statistical errors indicate one-sigma scatter in the six measurements per visit 
three images per exposure, two exposures). Compare to estimated errors in Table 

b Systematic uncertainty as tabulated in Table [| 



WF2 WF3 WF4 WF2 WF3 WF4 



Wavelength (A) 


Fig. 13.— Total background detected (DN/pixel in 1800 sec) in 
each of two exposures per filter obtained in the November (A) and 
December (o) 1995 visits. Bright galaxies (V 555 < 23 mag), stars, 
and bad pixels are excluded from this average. Filled symbols in¬ 
dicate the mean for each visit with error bars showing empirical la 
scatter around the mean, in good agreement with tabulated statis¬ 
tical errors (see Table pi). For this comparison, DN/pixel for each 
chip has been normalized to match the gain of the WF3 chip (see 
H95b). Decrease in the mean from November to December is in good 
agreement with the predicted change in the ZL flux due to viewing 
geometry. 


variation between the six measurements (three WF chips, 
two exposures per filter) from each visit. This scatter is 
well matched to the statistical errors that we predict based 
on our assessment of the errors accrued at each stage in the 
data reduction (see Table |3j) . Systematic errors are dom¬ 
inated by the flux calibration from DN to physical units 
and are not shown in the comparisons in Figure O as they 
are, of course, identical for all points. In Figure |F| we plot 
the total background flux measured in each bandpass in 


Fig. 14.— The total sky flux detected through each of the WFPC2 
filters in the November and December 1995 data sets, excluding 
bright galaxies (V 555 < 23ABmag). November 1995 results are indi¬ 
cated at the central wavelength of each filter; the December 1995 data 
are offset by +100A for clarity. The dashed line shows the ZL spec¬ 
trum at the appropriate flux; horizontal bars show ZL flux convolved 
with WFPC2 bandpass profiles. The total background detected with 
the FOS spectra are marked by o’s and offset left (November) and 
right (December) of the central F555W wavelength. Error bars on 
all points indicate one—sigma statistical errors (see Tables B| and p|). 
The floating error bars, arbitrarily plotted at 2500A, show system¬ 
atic uncertainties in the WFPC2 results at the indicated fluxes. The 
floating error bar at 7000A shows the systematic uncertainty in the 
FOS results. Solid curves show the effective F300W, F555W, and 
F814W bandpass profiles. 


the November and December visits in physical flux units 
(ergs s -1 cm -2 sr -1 A -1 ) with statistical and random la 
errors plotted separately. The measurements are summa¬ 
rized in Table Q. The zodiacal light spectrum is shown at 
the flux level we measure in our ground-based spectropho- 
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tometry (see §j|) to demonstrate that background observed 
from HST is clearly dominated by zodiacal light. 

6 . HST/FOS OBSERVATIONS AND DATA REDUCTION 

The FOS field of view during our parallel observations 
is determined by the WFPC2 coordinates and the roll an¬ 
gle of the telescope, which we carefully specified. Ground 
based images taken with the Swope 1 meter telescope at 
Las Campanas Observatory were used to assure that the 
FOS field of view was free of detectable sources brighter 
than V = 25 AB mag. 

Based on calculations using an early version of the FOS 
exposure-time calculator, the FOS/RED was used in par¬ 
allel observing mode and configured with the 1.0-PAIR 
aperture and the G570H disperser during our October 
visit, producing spectra with ~ 4.5A resolution from 4500 
to 6800A. No exposures were taken during transits of the 
South Atlantic Anomoly. These data were intended to 
provide a me asurement of the ZL by the method out¬ 
lined in § |2.2| . but the count rate from the sky in this 
configuration was comparable to the instrument dark rate 
(~ 0.01 DN/sec per diode). Accordingly, the FOS observa¬ 
tions during the November and December 1995 visits used 
the FOS/RED, A -1 aperture (3.66 x 3.71 arcsec 2 ), and 
G650L disperser, which provided spectra with ~ 12 diodes 
(~ 300A) per resolution element for an aperture-filling 
source (25A per resolution element for a point source) 
from 3800-7000A and a signal-to-noise ratio of almost 
20 for the night sky. While narrow band imaging as such 
is not useful for measuring the ZL flux by its spectral fea¬ 
tures, these data do provide better spectral resolution than 
WFPC2 images, and a second, independent measurement 
of the total background. We limit our discussion to the 
data from the November and December visits. One 1300 
sec FOS exposure was obtained per orbit during these vis¬ 
its, executing within the time-span of the 1800 sec WFPC2 
observations to avoid conflicts in writing the data to on¬ 
board recorders. 

Many of the pipeline calibration procedures (described 
thoroughly in the HST Data Handbook) contribute negli¬ 
gible errors for our purposes and were adopted directly. 
In order of their application, the pipeline procedures 
used include the conversion from counts to count-rate 
(cnt_CORR), flat fielding (flt_CORR), wavelength cal¬ 
ibration (wav_CORr) , and flux calibration (ais_CORR, 
APR_cORR, and tim_CORr) . The A 1 aperture defines the 
throughput standard for other apertures, so that the rela¬ 
tive correction (APR_CORR) for the A 1 aperture is unity. 
The G650L has no observed time variation, so that the 
correction factor applied by the TlM_CORR procedure is 
also unity. We were able to improve on the dark subtrac¬ 
tion procedure, and accordingly used our own methods for 
that step, as described below. Unfortunately, much of the 
information needed to obtain accurate surface photome¬ 
try with the FOS, namely the solid angle of the apertures 
and PSF of the instrument, has never been provided by 
STScI. We have therefore supplemented the pipeline cal¬ 
ibration with an explicit measurement of the PSF of the 
instrument. The final flux calibration has a systematic 
uncertainty of 4.5%. Errors at each step in the calibra¬ 
tion are summarized in Table |^, with the total systematic 
error tabulated as described in Below, we summarize 



pixel 

Fig. 15.— Average of 25 FOS dark exposures (1300 sec), demon¬ 
strating the stability of the pixel-to-pixel structure of the dark sig¬ 
nal. Strong features appear where diodes have been turned off due 
to poor performance. Sub-stepping along the diodes allows complete 
spectral coverage over these dead diodes. 


the calibration of the FOS data and briefly discuss the 
results. Greater detail can be found in Bernstein (1998). 
Once again, the following section is intended for readers 
who are interested in the details of the data analysis. 

6.1. Dark Subtraction 

The majority of the FOS instrumental background is 
Cerenkov radiation, caused by cosmic rays striking the 
photocathode, rather than thermal dark current. Conse¬ 
quently, the instrumental background can vary by factors 
of two between exposures, a fact ignored in the pipeline 
calibration. As with the WFPC2, the relative dark rate is 
a function of instrument geometry and is quite stable over 
the detector. The G650L low-resolution grating illumi¬ 
nates only 144 of the 512 diodes; therefore, the unillumi¬ 
nated portion can be used to identify the dark rate during 
a given exposure. On the suggestion of L. Petro (private 
communication), we produced a “superdark” using 25 dark 
frames taken in the HDF parallel program (Program 6342, 
Freedman; Program 6339, Petro) and scaled this to the 
level indicated by a subset of the unexposed diodes (pixels 
900-1100) in each of our spectra. The exposed and unex¬ 
posed regions of the diode array are labeled in Figure |l5], 
where we plot the averaged dark spectrum. To determine 
the accuracy of this dark-subtraction method, we have 
reduced darks which were taken during our own orbits ex¬ 
pressly for this purpose. The test reduction of these darks 
shows no systematic error in the dark subtraction by this 
method, and statistical errors are dominated by shot noise 
(±2.2 DN/pixel). The mean dark level can be determined 
in an individual frame to roughly ±0.3DN/pixel. With a 
sky signal of roughly 30 DN/pixel in the low-dispersion 
configuration, this dark-subtraction method introduces a 
0.5% random error to each spectrum overall and 0.4% sys¬ 
tematic uncertainty to each resolution element. 

6.2. Flux Calibration 

The surface brightness of an aperture-filling source is 
a function of the pipeline-calibrated spectrum, F( A) in 
ergs s -1 cm -2 A -1 , the detector solid angle, Cl, and the 
aperture correction, T(A-l)xP, which includes a term for 
flux lost at the A-l science aperture (T(A— 1)) and at the 
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Table 5 


FOS Background Flux: Errors per 
Resolution Element 


Random 

Systematic 

Poisson noise (§|fj) 

2 % 


Dark subtraction (§6 1) 

0.5% 

0.4% 

Fiducial standards (§ 4. 2.5|l 


1 % 

Point source flux cal. §312) 


1.5% 

Aperture correction (§3.2) 


2 % 

Solid angle (§|6.2|) 


2 % 

Cumulative** 

(±2.1%) 

[± 2.8%j 


Statistical errors have been combined in quadrature to ob- 
tain a cumulative, one-sigma error. Systematic errors have been 
combined assuming a flat probability distribution for each con¬ 
tributing source of error. The resulting systematic error is roughly 
Gaussian distributed, and the quoted value is the 68% confidence 
interval. For a detailed discussion see §0. 


detector ( D ): 


I(X) = 


F( A) T(A - 1) D 

n ' 


( 2 ) 


The pipeline calibration (calibration of the spectrum as 
appropriate for a point source) depends on the stability 
of the instrument for relative flux calibration and on the 
accuracy of the fiducial standards for absolute flux calibra¬ 
tion. The accuracy of the secondary standard star system 
does not dominate the FOS calibr ation uncertainty; it was 
already discussed briefly in §4.2.5 and is discussed further 
in Paper II. 

The pipeline calibration converts DN sec” 1 per diode to 
ergs sec” 1 A -1 , as appropriate for point sources, achiev¬ 
ing repeatability of 1-2% for point sources only if they 
are centered in the aperture to within ±0.2 arcsec (Keyes 
1995). This sensitivity to centering implies at least two 
serious complications for surface photometry. First, the 
transmission efficiency across the photocathode can vary 
by as much as 20% over surface areas corresponding to 
10 diodes. Second, and more importantly, the PSF is not 
well contained within either the aperture, the detector, or 
both. An accurate aperture correction is therefore crucial 
to surface photometry. 

The aperture dilution factor, T(ap) x D , given in the 
FOS Instrument Handbook is an estimate produced by 
modeling based on the OTA with post-COSTAR configu¬ 
ration. While no official Instrument Science Report exists, 
unofficial estimates for the monochromatic transmission 
of the A 1, post-COSTAR configuration at 6500A range 
from 97% (R. Bohlin, private communication) to 95% (The 
Data Handbook V3.0). Because this factor is crucial to our 
result, we have recalculated it using data taken for this 
purpose as part of the FOS calibration program (Proposal 
5262, Koratkar). These data were taken in ACQ/IMAGE 
mode, which uses no dispersing element, providing a two 
dimensional image in the diode plane. The stepping pat¬ 
tern used for the observations was kindly provided to us 
by E. Smith at STScI. Our reduction and analysis of these 
data are described in Bernstein (1998). We find that 98% 
of the flux from a point source is contained within the A 1 


aperture at the focal plane, and 96.5% of that flux is then 
imaged onto the 1.29 arcsec spatial extent of the diode 
array: T(A-1)=0.98 and D = 0.965. Thus, we find 

T(A -l)xfl + 0.002 ss 0.95, (3) 

in which we have included a small (0.002) correction for 
conversion from the “white” light of the ACQ/IMAGE 
data to our central wavelength of 5500A. The statistical 
error in this estimation is much less than 1%, but the 
systematic uncertainty may be as large as 2%. 

The A 1 science aperture measures 3.63 x 3.71 arcsec 2 , 
with 1% errors in both dimensions. In the spatial direc¬ 
tion the solid angle is determined by the diodes, which 
are 1.289(±1%) arcsec in the spatial direction. This val¬ 
ues is based on a laboratory measurement made before 
launch (Instrument Science Report ISR CAL/FOS-019), 
and has been corrected for the measured change in the FOS 
plate scale before and after COSTAR was installed (ISR 
CAL/FOS-123,141). The effective solid angle through the 
A-l aperture is then 3.63x 1.29 arcsec 2 , or 4.68 arcsec 2 with 
an uncertainty of ~ 2%. 

7. HST/FOS: RESULTS 

The combined averages of the six spectra taken dur¬ 
ing the November and December visits, respectively, are 
shown in Figure [h| With roughly 80 counts per diode, 
and twelve diodes per resolution element, the statistical 
error per resolution element for one spectrum is roughly 
4%. Statistical errors in the averaged spectra are roughly 
1.5-2%, indicated by the error bars which are placed one 
per resolution element (300A). The dark subtraction con¬ 
tributes an error of less than 0.5%. We do not show sys¬ 
tematic uncertainties, as they will affect both November 
and December data in the same way and are irrelevant 
for this comparison. We find that the integrated back¬ 
ground flux was brighter in November than in Decem¬ 
ber by 2 ± 0.5%, in good agreement with the expected 
change based on the increase in the path-length through 
the zodiacal plane between visits and empirical estimates 
from ground-based observations (Levasseur-Regourd & 
Dumont 1980). 


















16 


Detections of the Optical EBL: Results 



Fig. 16.— Surface brightness of the total background detected in 
the FOS spectra. The solid and dashed lines show the average of all 
six spectra taken during the November and December 1995 visits, 
respectively. The error bars indicate la statistical errors of 1.5-2% 
per 300A resolution element. 



Fig. 17.— Color of the detected background with respect to solar. 
Solid and dashed lines show the detected background in the Novem¬ 
ber and December 1995 data sets, respectively, divided by the Neckel 
& Labs (1984) solar spectrum at matched resolution. Straight lines 
show the linear fits. The error bars show the statistical error per 
resolution element in the FOS spectra. The adopted solar spectrum 
is simply the fiducial spectrum from which the ZL color is defined, 
and so it contributes no error here. For further discussion, see Paper 
II. 


Employing the usual defi nitio n of the ZL color relative 
to the solar spectrum (see §2.2), we find (7(7000,4000) = 
1.044 for the November visit, and 1.075 in the December 
data, as shown in Figure 17 with lcr statistical uncertainty 
of 0.05%. The reddening of the ZL relative to the solar 
spectrum is a function of the scattering angle and, thus, 
also of time of year. Again, the trend we observe is in 
good agreement with the Helios Space Probe observations 
and others (see Leinert etal. 1981, and Leinert etal. 1998 
for a review). We stress, however, that the FOS spec¬ 
trum includes both the EBL and the ZL, and therefore we 
cannot measure the color of the ZL explicitly from these 
observations. It is only possible to determine the ZL color 
separately from the EBL by explicit measurement of ZL 
absorption features over a wide range in wavelength. All 
published colors of the ZL to date rely on broad-band ob¬ 
servations. Separation of the ZL from other the EBL and 
DGL is discussed further in §11 and in detail in Paper II. 

An error budget for the absolute flux calibration of these 
data is shown in Table [| as the error per resolution ele¬ 
ment. The random error in the mean flux over the en¬ 
tire spectrum goes down as the square root of the number 
of resolution elements (x/TO). As discussed in § |ll] , ran¬ 
dom and systematic errors have been combined assuming 


Gaussian and flat distributions, respectively. For compari¬ 
son with the WFPC2 observations, we have integrated the 
FOS spectrum through the F555W bandpass (see Figure 


14) and find the two results in good agreement. This com¬ 


parison is subject to the systematic uncertainties in both 
data sets. 


8 . DIFFUSE GALACTIC LIGHT (DGL) 

8.1. Structured Component 

A diffuse, non-isotropic, optical background is produced 
in the Milky Way by scattering of the optical interstellar 
radiation field (ISRF) by Galactic dust. The same dust 
is heated by the UV ISRF, causing it to produce ther¬ 
mal IR emission. It is not surprising, therefore, that the 
thermal Galactic emission seen in the IRAS 100 /irn maps 
correlates well with the surface brightness of the optical 
diffuse Galactic light (DGL), as both are proportional to 
the column density of the dust and the intensity of the 
ambient ISRF along the line of sight. In Tables |(i| and [5J 
we give a representative summary of the observed corre¬ 
lations between optical and 100pm fluxes for regions with 
low to moderate 100 pm intensities (N(Hl)< 5x 10 20 cm -2 , 
Jioo < 5 MJy sr -1 ) and a range of Galactic orientations. 

As evident from those results, there is only moder¬ 
ate agreement concerning the exact scaling relations be¬ 
tween the optical DGL and thermal emission at any wave¬ 
length. Measurement errors in the IR, optical, and UV 
intensities are > 10% in most cases and are one cause 
for variations between results. However, asymmetry in 
the scattering phase function of Galactic dust also con¬ 
tributes to the variable scaling relations seen between dif¬ 
ferent lines of sight. Strong forward scattering causes 
lower optical surface brightnesses at both high latitudes 
(|6| > 50°) and at longitudes away from the Galactic Cen¬ 
ter (130° < l < 230°) (see Draine & Lee 1984, and refer¬ 
ences therein; Onaka & Kodaira 1991; Witt, Friedmann, 
& Sasseen 1997). Both trends are evident from the data 
shown in Tables || and [?]. 

While these results suggest a range of appropriate cor¬ 
relation factors, they do not identify a single appropriate 
scaling law for our purposes for two reasons. First, vari¬ 
ability in the measured IR-optical and UV-optical cor¬ 
relations is evident within a single cloud, as well as be¬ 
tween clouds (see, for example, Figures 5 & 6 in GT89 
and Table 1 in Laureijs, Mattila, & Schnur 1987). This 
suggests that the observed systems may be dense enough 
that self-shielding and further complications come into 
play. These clouds have been selected precisely because 
the optical and IR emission is bright enough to be readily 
observed: while the IR flux levels and N(Hi) column den¬ 
sities of the clouds listed in Tables |] and [] are low enough 
that the molecular gas fraction does not affect the cor¬ 
relation between dust column density (or extinction) and 
N(Hi), they are still roughly a factor of 10 higher than 
the values for our observed field, for which I wo ~0.4 MJy 
sr _1 (N(Hi) ~0.47xl0 20 cm -2 , or E(R — V)~ 0.009mag). 
Second, while empirical relations between the scattered 
and thermal DGL have have been published in the far- 
UV and at optical B- and R-bands, the expected surface 
brightness from scattering at 3000A is not clear from these 
results. Neither the optical depth of interstellar dust nor 
the ISRF is a monotonic function of wavelength between 
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Table 6 

Observed Correlations: /\(0.16/rm)//,,(!00/iin) 


Reference 

r a 

|5°| a 

I\ (0.16)//^(100) b 

7a(0.16)/C,(100) c 

G(100) d 

Witt et al. (1997) e 

(145) 

(30) 

72 

0 .86X10 -9 

2-8 

Hurwitz et al. (1991) f 

135 - 220 

> 40 

80(±10) 

0.96 x 10 -9 

1-5 

Sasseen etal. (1996) 

(270) 

(45) 

< 233 

< 2.8xl0 -9 

2-8 

Jacobsen et al. (1987) 

(70) 

(50) 

65(±25) 

0.78X10 -9 

1-2 

Witt et al. (1997) 

(290) 

(45) 

258 

3.10X10 -9 

2-8 


a Bracketed values of Galactic longitude (?) and latitude (b) indicate average co¬ 
ordinates for the data used in the calculation. 

h I\ is in photons s -1 cm -2 sr -1 A -1 . /„ is in units of MJysr -1 . 

C I\ is in ergs s -1 cm -2 sr -1 A -1 . 

d /„(100/im) is in MJysr -1 . 

e Witt et al. (1997) use a model which is based on their observations and which 
includes the average scattering angle and phase function along the line of sight to 
predict /A(0.16^m)//„(100/im) at the Galactic position indicated. 

f We calculate /A(0.16/im)//„(100/nn) for Hurwitz et al. (1991) from the points 
plotted in their Figures 2a and 2b. 


Table 7 

Observed Correlations: / z ,(A)//,,(100/iin) 


Reference 

in 

b(°) 

C(0.45)/C(100) a 

G(0.65)//„(100) a 

G(0.90)//„(100) a 

G(100) a 

GT89 - irl b 

174 

-42 

0.36X10 -3 

l.lxlO -3 


11.4 

Laureijs et al. (1987)° 

211 

-37 

0.48X10 -3 



6.3 

GT89 - ir2 b 

235 

37 

l.lxlO -3 

2 .2xl0 -3 

< 1.6xl0 -3 

3.6 

Paley et al. 1991 

104 

-32 

4.8xl0 -3 

8 .0xl0 -3 

11 .Ox 10 -3 

2.5 

GT89 - ir3 b 

38 

45 

2 .6xl0 -3 

4.4xl0 -3 

6 .Ox 10 -3 

5.9 


a /„ is in units of MJysr 1 , throughout. 1 MJysr 1 = 10 17 c/ A 2 ergs s 1 cm 2 sr 1 A A 
b Guhathakurta & Tyson (1989) - field designation. 
c From the data in Table 1 of Laureijs etal. (1987) 
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Table 8 

Model Correlations: /(A)//(100/im) 


A(jLzm) 

r A /N(Hi) a 

^ A b 

S(g, b) c 

7 A (A)//„(100) d 

I u ( A)/J„(100) e 

7 A (A)/G(100) f 

0.16 

0.120 

0.410 

0.42 

181.96 

1.86X10 -4 

2.18X10 -9 

0.25 

0.130 

0.550 

0.50 

94.54 

3.69X10 -4 

1.77X10 -9 

0.30 

0.105 

0.581 

0.54 

118.09 

7.97X10 -4 

2.66X10 -9 

0.45 

0.068 

0.600 

0.56 

68.75 

1.57x 10 -3 

2.32X10 -9 

0.55 

0.047 

0.600 

0.59 

42.85 

1.78X10 -3 

1.77X10 -9 

0.65 

0.042 

0.600 

0.61 

32.66 

2.24x10“ 3 

1.59X10 -9 

0.80 

0.032 

0.540 

0.68 

18.82 

2.41X10 -3 

1.13X10 -9 

0.90 

0.025 

0.500 

0.71 

11.76 

2.14X10 -3 

7.94xlO -10 


a Optical depth as a function of Hi column density in units of 10 2O cm -2 . Values are 
from Savage & Mathis (1979) and are in good agreement with the standard value of 
Ay/N(Hi)=0.06mag/(10 2 °cm -2 ) (Bohlin, Savage & Drake 1978). 

b Albedo values from models of Draine & Lee (1984). 

Scattering phase function in terms of asymmetry parameter, g and Galactic latitude (b = 50°), 
as calculated by Draine & Lee (1984). 

d / A (A) is in photons s - 1 cm - 2 sr - 1 A -1 . /j,(100/Ltm) is in MJysr -1 . 

e /„(A) is in MJysr -1 . J„(100^im) is in MJysr -1 . 

f / A (A) is in ergs s -1 cm -2 sr -1 A -1 . I„(100/im) is in MJysr -1 . 


1600A and 4500A (see Savage & Mathis 1979 and Mathis, 
Mezger, & Panagia 1983). 

To better understand the contribution of non-isotropic 
DGL over the full range of our observations, we have used 
a basic scattering model to predict the scattered light from 
dust. We then compare the results of this model to the 
observed DGL at UV and optical wavelengths. 

Assuming the Galactic cirrus along the line of sight in 
question is optically thin (extinction, A\ < 1.08mag), the 
surface brightness of scattered light off of interstellar dust 
can be expressed as 

h = [1 - l.li/Vsin b\. (4) 

in which j\ is the flux of the radiation field in ergs s -1 
cm -2 sr -1 A -1 ; lo\ is the effective albedo of the dust; r A 
is the optical depth; and the term in brackets is the back- 
scattered intensity in terms of Galactic latitude, b, and 
the average phase function of the dust, g (Jura 1979). For 
strong forward scattering, g ~ 1; for isotropic scattering, 
g ~ 0. We take the ISRF flux, j \, from the Mathis et al. 
(1983) estimate for the Solar Neighborhood (lOkpc from 
the Galactic center). As our observations are b = 60° 
from the Galactic plane and l =206?6 from the Galactic 
center, this estimate is probably slightly high. We take 
the dust albedo from the results of Draine & Lee (1984), 
which are based on an exponential distribution in grain 
sizes suggested by Mathis, Rurnpl & Nordsieck (1977). 

The optical depth of Galactic dust, r, is well known to 
correlate strongly with hydrogen column density (see Sav¬ 
age & Mathis 1979, Boulanger & Perault 1988, and refer¬ 
ences therein). It is not surprising, then, that the thermal 
emission, /ioo, also correlates well with hydrogen column 
density. While optical depth is a physical manifestation 


only of the column density of dust, /ioo is also affected by 
the strength of the ISRF. We therefore use the observed 
/ioo and Ji 0 o/N(Hi) as calibrated by Boulanger & Perault 
(1988) from the IRAS 100 /im maps to obtain an effec¬ 
tive optical depth for our observations as follows. Optical 
depth can be written as a function of optical extinction 
and dust column density as 


t\ = 0.921 


Rx 


N{m)/E(B-V) 


N( Hi), 


(5) 


in which R\ = A\/E(B — V) is the usual expression for the 
normalized extinction. Several groups find N(Hi)/E (B — V) 
between 48xl0 20 and 50xl0 20 cm 2 mag 1 from mea¬ 
surements of the Hi densities from 21 cm line emission 
strength and the reddening to globular clusters and star 
counts (Bohlin, Savage & Drake 1978, Burstein & Heiles 
1982, Knapp & Kerr 1974). To get an effective opti¬ 
cal depth (weighted by the ISRF held strength which 
is at issue for scattering), we use the relation found 
by Boulanger etal. (1996) for the low-column density 
regime (N(Hi) <5xl0 20 cm -2 ): /loo/N(Hi)cx 0.85 MJy 
sr -1 /(10 2 ° cm -2 )| The fluxes in our held are roughly 0.4 
MJy sr -1 , or 0.47x 10 2 ° cm -2 . The predicted scattered 
fluxes from this model are shown in Table |[ Scattering 
angle is not considered in this model. Consequently, this 
estimate is conservative in the sense that it should over¬ 
predict the DGL for our observations, as the line of sight 
to our held is away from the Galactic center and the dust 
is forward scattering. 

3 A slightly different scaling, boo/N (Hi) oc 0.53 MJy 
sr —1 /(10 2 ° cm -2 ), is seen in from the DIRBE results (Boulanger 
etal. 1996). The difference is attributed to a well known calibration 
offset in the IRAS maps. Since we are using IRAS fluxes, we use the 
IRAS correlation. 
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This scattering model reproduces the observed flux ra¬ 
tios with reasonable accuracy in the range 1600-4500A 
(see Tables [?] and |(] at b > 45). The phase function changes 
by less than 10% at latitudes | 6 | > 50°, so the values shown 
in Table | 8 [ for which \b\ = 50° was used, are generally rep¬ 
resentative for high latitude fields. However, as noted by 
GT89, optical colors {B—R) and ( R—I ) are redder than 
a basic scattering model predicts. GT89 find values of 
I V {R)/R{B) = 3.2, 2, and 1.7 and I V {I)/I V {R) = 2.3, 2.1, 
and < 1.5 in three different fields. By comparison, the 
ratios we predicted are I V {R)/I V (B)=\A and I U (I)/R(R) 
= 0.95. A significant Ha contribution as the explanation 
for the red colors is ruled out by GT89. Variable scat¬ 
tering asymmetry with wavelength is another possible ex¬ 
planation, but strong wavelength dependence in the range 
4500-9000A has never been observed in the lab or in space 
(Witt etal. 1997, Onaka & Kodaira 1991, Laureijs etal. 
1987). The most plausible explanation is suggested by 
observations of reflection nebulae, which have high N(Hl) 
and show red fluorescence from molecular hydrogen, hy¬ 
drogenated amorphous hydrocarbons, and polycyclic aro¬ 
matic hydrocarbons. The relevance of such contributions 
to fields with 10 times lower N(Hi) and Iioo, as is the 
case for our data, is not clear, as the density of molecu¬ 
lar gas correlates only with high column densities, N(Hl) 
> 5xl0 2 O cm -2 . The results of GT89, in fact, do show 
that the degree of reddening is well correlated to the av¬ 
erage /ioo emission but not structure within the cloud. 
Self-shielding, local optical depth and local ISRF may be 
responsible for strong variations in the correlation between 
color and molecular gas density both in and between fields 
(Stark 1992, 1995). It seems conservative, therefore, to 
adopt optical colors found for the fields with the lowest IR 
flux in the GT89 sample, listed in Table [?]. Note that the 
IR flux in the 2 lower flux fields (denoted “ir2” and “ir3”) 
is still more than a factor of 10 higher than in our own. 

In summary, we estimate the optical flux in our field 
using our scattering model for A < 4500A, and adjust the 
predicted scattering model at redder wavelengths to match 
the average colors observed by GT89: I V (R)/R(B) ~ 1.8 
and I V {I) / I V {R) ~ 2.0. We apply this correction in the 
sense of increasing the long wavelength fluxes over that 
predicted by our models, so that the DGL estimate we 
use is, if anything, higher than is appropriate, although 
given the small total flux associated with the DGL even 
a large fractional decrease in our estimate of the DGL 
would have a negligible impact on our EBL results. The 
resulting spectrum is flat in I\, with a value of roughly 
0.9-1.0xl0~ 9 ergs s -1 cm -2 sr -1 A -1 from 3000-9000A. 
We note, also, that our scattering model was not depen¬ 
dent on Galactic longitude, which, again, makes ours a 
conservative overestimate of the DGL contribution to the 
total sky background, and our measurement of the EBL, 
therefore, a conservative underestimate in this regard. 

8.2. Isotropic Component 

Line emission and continuum processes from warm ion¬ 
ized gas in the Galaxy also contribute an isotropic com¬ 
ponent to the DGL. For \b\ > 5°, Reynolds (1992) finds 
that Ha emission strength matches the prediction of a 
path-length through a slab model for the galaxy, /(Ha) ss 
2.9 x ICR 7 esc |5| ergs s _1 cm -2 sr -1 . Fortunately, Ha 
emission, specifically, is irrelevant for us because the rel¬ 


ative throughput of the F555W bandpasses at Ha (~ 
6562A) is only ~ 10 % of the peak filter throughput. The 
strongest Ha emission expected in our field would con¬ 
tribute 0 . 01 xl 0 ~ 9 ergs s _1 cm -2 sr -1 A -1 , which corre¬ 
sponds to 0 . 01 % of the total background, and 1 % of the 
expected EBL. The next strongest line, [Om] at 5007A, is 
near the peak of the F555W sensitivity, but it is fainter 
than Ha by a factor of 20 (Reynolds 1985, Shields etal. 
1981) and will contribute at most 0.05% of the expected 
EBL. 

More important than line emission is the two-photon, 
free-free, and bound-free continua emitted by ionized gas 
with the density implied by the detected Ha emission. 
The combined spectrum of free-free, bound-free, and two- 
photon emission was calculated by Aller (1987) as a func¬ 
tion of electron and ion densities, and has been expressed 
by Reynolds (1992) as a function of the observed Ha emis¬ 
sion; it is a function of the temperature of the warm 
ionized medium. For our purposes, a conservative esti¬ 
mates of the isotropic continuum from gas with tempera¬ 
ture T ~ 10 4 is given by Aller (1987) and Reynolds (1992) 
as I X (< 3700A) < 0.3xl0 _9 ergs s _1 cm -2 sr -1 A -1 and 
I\(> 3700A) < 0 . 01 xl 0 ~ 9 ergs s _1 cm -2 sr -1 A -1 (see 
Aller 1987 and Reynolds 1992 for discussion). We include 
these contributions in our estimate of the DGL. 

9. MEASUREMENT OF ZODIACAL LIGHT 

To measure the ZL contribution to the background flux 
that we measured with HST, we obtained ground based 
spectra at 3900-5100A within the field of view of our 
WFPC2 images using the Boiler & Chivens spectrograph 
on the duPont 2.5m Telescope at Las Campanas Obser¬ 
vatory in Chile on the nights of 26-29 November 1995, 
concurrently with the HST observations of that field on 
29 November 1995. We have used those spectra to mea¬ 
sure the absolute flux of the ZL at 4650A with a precision 
of 0 . 6 %, and a systematic uncertainty of 1 . 1 %, using the 
method outlined in §2. We also mesure the color of the 
zodiacal light to be (7(5100,3900) = 1.05 ± 0.01. That 
measurement is discussed in detail in Paper II. 

To identify the ZL flux contributing to the FOS and 
WFPC2 measurements of the total sky flux, we need an 
absolute, flux-calibrated spectrum of the ZL from roughly 
2000-1/xm (see WFPC2 bandpasses plotted in Figure |14|). 
We obtain this spectrum by scaling a solar spectrum to the 
surface brightness value we measure for the ZL at 4650A 
and applying a small reddening correction redward and 
blueward of 4650A to compensate for the changing scat- 
tering efficiency of the interplanetary dust (see § | 2 . 2 | for 
discussion). Although our LCO measurement of the ZL 
color is quite accurate, it only covers a small fraction of 
the total wavelength range covered by our WFPC2 obser¬ 
vations. We therefore use our FOS observations to iden¬ 
tify the appropriate reddening correction, which cover the 
wavelength range 4000-7000A. Published measurements of 
the ZL color have absolute uncertainties as large as 10% 
and are in poor agreement with each other (see Leinert 
1998). Our FOS observations, by contrast, are accurate 
to better than 2 % in relative flux as a function of wave¬ 
length and are identical in line of sight and epoch to the 
WFPC2 and LCO observations. Like most measurements 
of the ZL, however, the FOS observations clearly include 
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R(mask)/r iso 


Fig. 18. — The ensemble flux recovered from apertures around 
detected galaxies (V > 23 AB mag) as a function of the radius of the 
apertures. Radius is given in terms of the limiting isophotal radius 
(’"iso) for each galaxy. The solid line shows the flux recovered when 
ri S o is unconstrained; for the faintest objects, r*i so may be less than 
3 pixels (the 90% encircled energy radius of WFPC2). The dotted 
line shows the flux recovered when a 3 pixel minimum is imposed for 
riso- Note that the curve has not yet converged at the limit of the 
mask size plotted here. 


the EBL flux. This fact prevents the FOS and LCO re¬ 
sults from unambiguously determining the color of the ZL. 
Fortunately, this is not an insurmountable problem: the 
WFPC2 observations can be used as a further constraint, 
as we describe in the following section. 

10. EBL FROM RESOLVED SOURCES IN WFPC2 IMAGES 
10.1. Ensemble Photometry 

The total flux from resolved sources defines a lower limit 
to the EBL. Typically, such minima are obtained by mea¬ 
suring the flux in individual resolved sources using stan¬ 
dard photometry packages, such as FOCAS (Valdes 1982, 
Jarvis & Tyson 1981) or SExtractor (Bertin & Arnouts 
1996), and summing the flux in the resulting catalog of 
objects (c.f. Tyson etal. 1989, Madau etal. 1996, 1998). 
While it is quite straightforward to measure an isophotal 
magnitude, it is very difficult, if not impossible, to mea¬ 
sure the total flux of an object including that flux in re¬ 
gions fainter than the noise level of the local sky. Efforts at 
measuring galaxy counts often attempt some sort of “to¬ 
tal magnitude” correction, either by scaling the isophotal 
aperture to include roughly twice the isophotal extent (see 
Bertin et al. 1998 for discussion) or by applying a magni¬ 
tude correction to galaxies near the detection limit (c.f. 
Srnail etal. 1995). 

Fortunately, individual source photometry is not nec¬ 
essary for estimating a lower limit to the EBL in our 
data. We have developed a simplified method of aper¬ 
ture photometry with which we can measure the flux from 
the ensemble galaxy population as a whole. We first use 
SExtractor to identify detectable sources and their isopho¬ 
tal radii (r; so ) using detection and extraction parameters 
very sim ilar to those employed in Williams etal. (1996) 
(see §4.4). We can then measure the “sky” surface bright- 
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Z-^core (AB mag arcsec -2 ) 

Fig. 19.— Isophotal surface brightness, /Xi so , versus core surface 
brightness (mean surface brightness in the brightest 4 pixels), fi c ore, 
for all detected galaxies in the EBL images in unit magnitude bins 
22 < V 555 < 28; the faintest bin with data shows 27 < V < 28. 
In each plot, the diagonal lines (right) mark the limit /q so = //core; 
the horizontal lines (left) show the mean /q so , defined by sky noise 
(lcr s ky)- As expected, this does not change with magnitude. The 
distribution in /x c ore narrows at fainter apparent magnitudes, sug¬ 
gesting that surface brightness biases and detection limits may be 
causing incompleteness at the fainter magnitudes. The detection 
limit occurs at V = 27.5 AB mag, at which the average galaxy has 
Miso — A^core ~ 1-5 mag arcsec -2 . The requirement for detection in 
our data is fi C ore < 2cr s k y in the central 4 pixels. 


ness, /i s ky, for an image by masking all detected sources 
(V 555 < 27.5 AB mag) and computing the average surface 
brightness of the remaining pixels. In doing so, we im¬ 
plicitly assume that extragalactic sources fainter than this 
limit contribute negligibly to the EBL, so that the appar¬ 
ent background level in the image consists only of fore¬ 
ground sources f| To estimate the total surface brightness 
from extragalactic sources and foregrounds, p. 0 bj+sky, we 
mask all stars and only galaxies brighter than I4 u t = 23AB 
mag and again compute the average surface brightness of 
the remaining pixels. We then isolate the flux from re¬ 
solved sources, the minimum value of EBL23, by differ¬ 
encing the two surface brightness estimates: minEBL23 
= p-obj = Mobj+sky — /xsky- In contrast to standard aper¬ 
ture photometry, /x 0 bj+sky will always include all of the 
light from faint galaxies. However /x s k y can include galaxy 
flux if the masks are too small. 

Not surprisingly, minEBL23 from this method of en¬ 
semble photometry is a strong function of the mask size. 
The growth curve plotted in Figure ^ 8 | shows that at least 
20% of the flux from galaxies within 4.5 magnitudes of 
the detection limit (23 < V < 27.5 AB mag) lies beyond 
t/2 x riso for individual galaxies. As discussed in Williams 
etal. (1996), it is possible for the faintest galaxies detected 
to have isophotal radii which are smaller than the WFPC2 

4 We point this out not because we believe that the flux in galaxies 
beyond the detection limit is negligible, but rather as a reminder that 
the sum in detected galaxies is by definition a minimum value for the 
detected EBL. 
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Fig. 20.— The same quantities as plotted in Figure |19| are shown 
for the galaxies in the HDF F606W catalog (Williams etal. 1996). 
The detection limit of the HDF catalog for Vqqq is quoted as 28.2 
AB mag, which corresponds to a 10cr detection in the drizzled images 
2cr s ] cy in 4 pixels in the original images). At V 606 = 28.2AB mag, 
a typical galaxy has (^i so — Mcore) < 1 mag arcsec -2 . See Appendix 
B.2 for discussion. 


90% encircled energy radius for a point source (0.3 arcsec, 
or 3 pixels). To explore the effect of this resolution limit, 
we also show in Figure |B the growth curve which results 
if we impose 0.3 arcsec as a minimum for r- lso and increase 
the mask size relative to that new starting radius. Natu¬ 
rally, the improvement in recovered flux is significant: in 
the HDF catalog of Williams etal. , nearly all galaxies 
with Veoe > 27.5 ABmag are affected by this size require¬ 
ment. While the instrumental PSF clearly contributes to 
the increasing recovered flux with radius, the growth curve 
is not entirely due to PSF effects: the 100% encircled en¬ 
ergy radius would be reached by a mask of radius 2.5ri so if 
the smallest galaxies were physically contained within the 
original isophotal radii. 

In addition to the fact that apertures extending to twice 
each galaxy’s isophotal area will not recover all of the 
galaxy light from detected galaxies, the “local sky” es¬ 
timates used in standard photometry packages come from 
regions just beyond the apertures for each source. These 
sky estimates will undoubtedly include a significant frac¬ 
tion of this missed light, doubly compounding the pho¬ 
tometry errors. 

While the growth curve in Figure |l8| is beginning to 
level off by 4ri so , it is clearly not flat. The covering factor 
of galaxies to 4r; so is roughly 20% in these images, and 
roughly 30% of the galaxies defined by these apertures 
overlap. In the HDF images, more than 1800 galaxies are 
found beyond our detection limit, with a covering factor 
of roughly 80% to 4rj so . These facts suggest that we have 
reached a confusion limit of sorts, in the sense that the 
wings of detectable galaxies are overlapping. Moreover, 
these facts suggest that the wings of detectable galaxies 
contribute flux to a significant fraction of the pixels used 
to measure the foreground “sky” level, producing an ex- 


tragalactic “pedestal” level which blends with the diffuse 
flux from ZL and DGL. Even though the flux level of our 
growth curve in Figure E is clearly converging, it cannot 
be used to quantify this pedestal of overlapping galaxy 
wings. 

To estimate the flux contained in the wings of galaxies 
beyond 4rj so in our own images, we have constructed a 
Monte Carlo simulation of the contribution to a random 
point on the sky from randomly distributed, detectable 
galaxies. The surface number density of detectable galax¬ 
ies was drawn from the published HDF catalog and ex¬ 
ponential light profiles were adopted for all galaxies, with 
scale lengths and central surface brightness matching the 
galaxies in the EBL and HDF images. Efforts to char¬ 
acterize the light profiles of faint galaxies have produced 
evidence for both exponential disks and flatter, irregular 
profiles (for example, see Smail etal. 1995, Driver etal. 
1995, Brinchmann etal. 1998, Driver etal. 1998 and ref¬ 
erences therein). For the faintest galaxies, profiles are un¬ 
constrained. We adopted exponential profiles here, in part 
because they produce the most conservative (smallest) es¬ 
timate of the light beyond the 4ri so apertures. Also, given 
that the physical scales at those radii are large, we expect 
to be beyond any central bulges. We also note that the 
additional flux identified by our models from the faintest 
galaxies is less than 10% of the flux in galaxy counts by 
standard methods, and so is the adopted profiles for the 
faintest galaxies are not a critical issue to the minimum 
EBL estimates. The simulations are described in Ap¬ 
pendix B. For the F555W EBL images, we estimate that 
the additional flux from galaxies beyond 4ri so is roughly 
l.lxlO -10 ergs s -1 cm -2 sr -1 A -1 . The simulations are 
affected only by the surface brightness limits, galaxy pro¬ 
files, and the galaxy surface densities adopted. Similar 
levels are therefore found for the F814W images, for which 
galaxy parameters and detection limits are very similar to 
those at F555W. The F555W aperture masks have been 
used to recover the flux from sources in all three bands for 
reasons discussed below. In the very low signal-to-noise 
ratio F300W images, only 20% of the F555W sources are 
detected, and the F555W aperture masks for those which 
are detected extend to many more than four times the 
F300W isophotal detection areas. Due to the larger statis¬ 
tical uncertainties in the flux recovered by ensemble pho¬ 
tometry at F300W (see Table [Toj), simulations of the flux 
beyond the detection apertures were not warranted. 

By using different values for Wut (23, 24, 25 and 26 AB 
mag), we can isolate the flux contributed from galaxies 
in successive magnitude bins. Comparing these measure¬ 
ments with the int egrated flux from standard photometry 
methods (see §L4), we find that roughly 15%, 25%, 45% 
and 65% of the total flux in successive 1 magnitude bins 
between 23 < V555 < 27 is contained at radii between 
\/2riso and 4ri so . In our data, galaxies with magnitudes 
23 < V555 < 27.5 have 1.5 < A/i < 3.5 mag arcsec -2 
as shown in Figure |9[ These measurements are in broad 
agreement with the models in the literature for the frac¬ 
tion of flux which can be recovered as a function of A/i 
based on extrapolation of simple exponential profiles (c.f. 
Disney & Phillipps 1983, Davis 1990, Dalcanton 1998). 

The total flux in resolved galaxies as measured by en¬ 
semble photometry has two obvious implications for galaxy 
counts. First, the total flux in galaxy counts based on stan- 
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Table 9 

Summary of Measurements 


Source 

Bandpass 

Data source 

Flux 

Random 

Systematic 

Total 

F300W 

hst/wfpc2 

33.5 

(± 4.9%) 

[± 5.6%] 

Background 

F555W 

hst/wfpc2 

105.7 

(± 0.3%) 

[± 1.4%] 


F814W 

hst/wfpc2 

72.4 

(± 0.2%) 

[± 1.496] 


F555W a 

hst/fos 

111.5 

(± 0.7%) 

[± 2.8%] 

Zodiacal 

4600-4700A 

LCO 

109.4 

(± 0.6%) 

[± 1.196] 

Light 

F300W 

LCO b 

28.5 

(± 0.6%) 

[-1.1%,+1.2%] 


F555W 

LCO b 

102.2 

(± 0.6%) 

[-1.1%,+1.1%] 


F814W 

LCO b 

69.4 

(± 0.6%) 

[-1.3%,+1.1%] 

Diffuse 

F300W 

DGL model 

1.0 


[+2596,-50%] 

Galactic 

F555W 

DGL model 

0.8 


[+2596,-50%] 

Light 

F814W 

DGL model 

0.8 


[+2596,-50%] 


Note. —All fluxes are in units of 1x10 9 ergs s 1 cm 2 sr 1 A 1 . For a source with constant 
flux in F\, filters F300W, F555W, and F814W have effective wavelengths Ao(AA) =3000(700), 
5500(1200), and 8100(1500)A. For a source with a solar spectrum, effective wavelengths are Ao = 
3200, 5500, 8000A. 

a Observed FOS spectrum, convolved with the WFPC2/F555W bandpass to allow direct com¬ 
parison with the WFPC2 results. 

b LCO measurement of zodiacal light, extrapolated to the WFPC2 bandpass by applying a 
correction for changing zodiacal light color with wavelength relative to the solar spectrum. The 
zodiacal light flux through the WFPC2 bandpasses was identified using SYNPHOT models, the 
uncertainty due to which is included in the uncertainty for the filter calibration and is shared with 
the systematic uncertainty for the total background flux. 


Table 10 

EBL Results and Uncertainties 


Bandpass 

Random 

Systematic 

Combined 

EBL(±lcr) 


CTR (68%) 

eg (68%) 

a (68%) 



Detected EBL23 (WFPC2 

+ LCO) a 




F300W 

2.1 

1.5 

2.5 

4.0 (±2.5) 

F555W 

0.6 

1.3 

1.4 

2.7 (±1.4) 

F814W 

0.4 

0.9 

0.0 

2.2 (±1.0) 

Minimum EBL (WFPC2) a 




F300W 

0.19 

0.13 

0.22 

3.2 (±0.22) 

F555W 

0.003 

0.009 

0.01 

0.89 (±0.01) 

F814W 

0.002 

0.007 

0.01 

0.76 (±0.01) 

Detected EBL23 (FOS ± LCO) a 



F555W 

0.7 

2.7 

2.8 

8.5 (±5.6) 

Flux from detected sources 

in HDF (m 

> 23 AB mag) 


F300W 

F450W 

F606W 

F814W 




0.66 

0.51 

0.40 

0.27 


Published number counts 13 
F300W (18 < U300 < 23 AB mag) 
F555W (15 < E 555 < 23 AB mag) 
F814W (13 < /si 4 < 23 AB mag) 


0.27 (±0.05) 
0.49 (±0.10) 
0.65 (±0.13) 


Note. —All fluxes and errors are given in units of 10 9 ergs s 1 cm 2 sr 1 A 1 . 

a The systematic and statistical errors have be.en combined assuming a flat or Gaussian proba¬ 
bility distribution, respectively, as discussed in § |ll[ We equate ler combined errors with the 68% 
confidence interval, as the combined errors are nearly Gaussian distributed. Individual sources of 
error contributing to these totals are summarized in Tables 3 and 4 of this paper and in Table 1 
of Paper II. 

b Estimated errors correspond to uncertainties in the fits to published galaxy counts. The values 
given correspond to 0.081 x 10 _2 °, 0.46 x 10 _2 °, and 1.5 x 10 _2 ° in units of ergs s -1 cm -2 sr -1 
Hz -1 and are consistent with those used in Pozzetti et al. (1998). 
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dard photometry will significantly underestimate the EBL. 
This point is discussed further in the following sections. 
Second, because standard photometry will miss fraction¬ 
ally more light from the faintest galaxies, the galaxy counts 
which result will have an artificially shallow slope at the 
faint end. Our results allow us to derive “aperture correc¬ 
tions” for faint galaxy photometry as a function of central 
and isophotal surface brightness. This is done in Paper 
III. We can use these corrections to recalculate the surface 
number density of galaxies as a function of magnitude. 
As we discuss in detail in Paper III, the corrected number 
counts do not flatten out at the faintest limits of the HDF. 


10.2. Discussion: Minimum EBL23 


Using the ensemble photometry method described above, 
we measure the flux from detectable sources relative to the 
mean sky level in the F300W, F555W, and F814W im¬ 
ages of our own data to be 3.19xl0 -9 , 6.02xl0~ 10 , and 
5.16xl0 _1 °, respectively^] For the F300W and F814W 
images, we use the masks derived from the F555W images 
in order to guarantee that the same sources are contribut¬ 
ing to the minimum EBL23 at all wavelengths. The 1800 
galaxies detected in the HDF with Vqoq magnitudes in the 
range 27.5-30 AB mag (see Figure ||) are clearly not in¬ 
cluded in our estimate of the minimum EBL23 derived by 
ensemble photometry. 

The contribution from galaxies in the HDF catalog with 
V 555 > 27.5 AB mag is 0.57x 10 _lo ergs s -1 cm -2 sr -1 A -1 
at V 555 as measured by standard photometric methods. In 
keeping with the discussion in the previous section, we es¬ 
timate that only 35% of the light is recovered from galax¬ 
ies with 27.7 < Vgo 6 < 30ABmag in the HDF, so that 
the true flux from these sources is roughly 1 . 8 x 10 _1 °ergs 


,-1 


s~" cm -2 sr -1 A -1 . Adding this to the flux from de¬ 
tected sources (23 < V 555 < 27.5) in our EBL field, we 
find the total flux from detected sources (23 < V 555 < 30) 
to 4r iso apertures to be 7.8x 10 _1 °ergs s -1 cm -2 sr -1 A -1 . 
Adding to this the estimated extragalactic contribution to 
the background sky level beyond 4ri so (Figure ]B2|), we 
identify the minimum flux from detected galaxies with 
U 555 > 23ABmag to be 8.9xl0 _lo ergs s -1 cm -2 sr -1 
A" 1 . We emphasize that this estimate of the minimum 
EBL23 is indeed a minimum estimate from which sources 
will be excluded due to the surface brightness biases and 
detection limits of our own images as well as the HDF (see 
Figures ^9] and ^0|). Incompleteness due to surface bright¬ 
ness detection limits is discussed further in Paper III. 

Following the same method for F300W and F814W data 
leads to the minimum EBL23 values summarized in Table 
10, The total combined statistical and systematic error for 
this minimum EBL23 measurement is roughly ±lxl0 -11 
ergs s -1 cm -2 sr -1 A -1 , xlOO smaller than the error for 
the EBL23 detections because the conversion to physical 
units can occur after foregrounds are subtracted. Com¬ 
paring the total flux from detected galaxies measured us¬ 
ing ensemble photometry versus standard methods, the 
photometry errors affecting standard methods are clearly 
worse at UV wavelengths where signal-to-noise ratio is 
lower, calibration is less accurate, and the instrumental 

5 For comparison, the flux measured by standard photometry for 
sources with AB magnitudes between 23 and 27.5 in the HDF are 
6.7, 5.1, 3.8, and 2.5x 10 _10 ergs s -1 cm -2 sr -1 A -1 in the F300W, 
F450W, F606W, and F814W, respectively. 
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Fig. 21.— Filled circles show the EBL23 obtained using the 
HST/WFPC2 EBL images, the ZL flux, measured at (4650A), and 
the nominal ZL color as described in §lll The solid error bars (cen¬ 
tered on the filled circles) show the combined lcr errors while the dot¬ 
ted (offset by +100A) and long-dashed (offset by -lOOA) error bars 
show lcr systematic and random errors, respectively. The hatch- 
marked region shows the lcr uncertainty in the detected EBL due 
to uncertainty in ZL color. The lower limit arrows connected by a 
dashed line indicates the total flux from individually photometered 
galaxies with magnitudes 23 < V 555 < 30 AB mag in the HDF cata¬ 
log. The u-shaped lower limit arrows show minEBL23, the flux from 
detectable galaxies in the EBL fields with 23 < V 555 <_28ABmag 
as determined by ensemble photometry as described in §| 10 [ 


PSF is broader than in the other filters. 


11. EBL DETECTIONS 


The detected EBL in each of the three WFPC2 bands 
is the difference between the total sky flux in the WFPC2 
images, /Total, a nd t he foreground flux from ZL and DGL. 
As described in § 3T, the DGL contribution has no spectral 
features, is relatively flat in I\, and is easily subtracted. 
The spectrum of the zodiacal light is uniquely defined by 
our measurement of the mean flux of the ZL at 4650A, the 
solar spectrum, and the color of the ZL as follows: 


/zl(A)=J q (A) <7(A, 4650) 


/zl(4650) 
J 0 (465O) ' 


( 6 ) 


(See the Appendix A for a discussion of the solar reference 
spectrum, used and the convolution of /zl (A) with the 
WFPC2 bandpasses to obtain the ZL flux through each 
filter.) 

Both the mean flux and color of the zodiacal light in the 
range 3900-5100A are uniquely determined by our LCO 
data. However, the color over the wavelength range of 
the WFPC2 data is less certain. From our LCO data, 
we find (7(5100,3900) = 1.05 ± 0.01 from absolute sur¬ 
face spectrophotometry of the ZL. As discussed in §[?] we 
find (7(7000,4000) = 1.044 ± 0.01 from our FOS spectra. 
Although this is in excellent agreement with the LCO re¬ 
sults, extragalactic contributions are included in the FOS 
spectrum. We have therefore devised a method to further 
constrain the ZL color using the WFPC2 results. 
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Uncertainty (xl0~ 9 ergs s~‘ cm -8 sr -1 A -1 ) 
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Uncertainty (xlO -9 ergs s _1 cm' 2 sr -1 A -1 ) 


Fig. 22.— Combined probability distributions of the errors (sys¬ 
tematic, random, and combined) for the EBL flux detected in the 
F300W (panel a), F555W (panel b), and F814W (panel c) WFPC2 
bandpasses. A flat (Gaussian) probability distribution was assumed 
for each contributing systematic (random) error. The vertical lines 
show 68.3%, 95.4%, and 99.7% confidence intervals, which we use to 
define 1, 2 and 3cr values for the combined errors. 


Between 2000A and 1 /jm. the color of the zodiacal light 
relative to the solar spectrum as a function of wavelength 
is empirically known to be roughly linear with wavelength; 
the mean flux with respect to the Sun becomes redder by 
roughly 5%/1000A over the range from 2500A-1.4/rm. No 
deviations from a linear color are apparent in our FOS 
spectra, which are the first observations to have the re¬ 
quired broad-band calibration sensitivity to address this 
question. To constrain the ZL and EBL surface bright¬ 
ness in all three bands, we therefore begin by assuming 
that the color of the ZL is a linear function of the solar 
spectrum with wavelength and that any small deviations 
from linearity average out over the 1000A bandwidths of 
the WFPC2 filters. 

Only if the EBL is flat in I\ would the color of the 
EBL plus ZL be the appropriate color to adopt for the 
ZL. If, for example, the EBL is blue in I \, then the ZL 
is slightly redder. At 3000, 5500, and 8000A, we have 
already identified minima for the EBL. If we adopt these 
minima for the EBL at 3000 and 8000A, we can obviously 
infer that the maximum possible flux for the ZL is what 
remains when we subtract the DGL and the minimum EBL 
from the total surface brightness detected in each band. 
By considering the maximum possible Jzl at 3000 and 



Fig. 23.— The detected background light from galaxies fainter 
than V 555 = 23 AB mag is marked by large filled circles with 2cr er¬ 
ror bars. The error bars are dotted where they extend below the flux 
recovered from galaxies in the rang£_23 < V 555 < 28 AB mag by the 
aperture photometry method (see §flO|) shown by the u-shaped lower 
limit symbols. The integrated flux from galaxies in the HDF, from 
23 AB mag to the detection limit at each wavelength (29-30 AB mag), 
is marked by lower limit arrows connected by the thick, dashed line. 
The total flux in galaxy counts, HDF counts plus ground-based 
counts brighter than V 555 = 23ABmag is shown by lower limit ar¬ 
rows connected by a thin dashed line. The open circles (displaced 
by 200 A for display purposes) show the EBL we detect, plus the in¬ 
tegrated ground-based counts brighter than V 555 = 23 AB mag. We 
have not included uncertainties in the ground-based counts in the 
error bars shown. 


8000A, and comparing those maxima to the absolute flux 
which we have measured at 4650A, we constrain the color 
of the ZL to be 3.0%/1000A < C( A,A 0 ) <5.0%/1000A. 
This estimate is in excellent agreement with the color we 
find for the total background as measured by FOS, which 
was 4.4%/1000A redder than solar. Using these extrema 
for the color of the ZL, we can constrain the EBL at all 
three WFPC2 wavelengths (see Figure ^l|). We adopt this 
range as the 3er uncertainty in the ZL color. 

We calculate final statistical and systematic uncertain¬ 
ties from the WFPC2 and LCO measurements assuming 
a Gaussian probability distribution for each source of ran¬ 
dom error listed in Tables |] and |] and in Table 1 of Paper 
II. The final systematic uncertainties are calculated assum¬ 
ing a flat probability distribution for the systematic un¬ 
certainties listed in the same tables. Cumulative random 
and systematic errors are quoted and plotted separately to 
give a sense of the measurement limitations. In general, 
random errors reflect limitations in instrument sensitiv¬ 
ity and stability, while systematic errors are dominated by 
stability of the flux calibration of each instrument and the 
accuracy of the flux calibration for the instruments rela¬ 
tive to each other. Uncertainty in the DGL subtraction 
is not explicitly included in these combined errors, as the 
total DGL flux is relatively small and the error would not 
contribute significantly to the total error. The DGL esti¬ 
mate was intended to be conservatively large, so that the 
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EBL is, if anything, underestimated due to errors in the 
DGL. 

We have also calculated combined errors (systematic 
and statistical) to obtain a final confidence interval. For 
this we have again assuming a Gaussian probability dis¬ 
tribution for all of the random sources of error and a flat 
probability distribution for the systematic uncertainties. 
The combined error probability distributions are shown in 
Figure [2^. In general, non-Gaussian systematic errors and 
Gaussian random errors cannot meaningfully be combined 
in quadrature. Nonetheless, given the large number (16) 
of individual sources of uncertainty contributing to our 
final errors, the combined uncertainty has a nearly Gaus¬ 
sian distribution, as expected. We therefore equate the 
lcr combined errors for each bandpass with the 68 % con¬ 
fidence intervals for the mean EBL23 detections and the 
lower limits, minEBL23. The final errors are summarized 
in see Table [io]. 

12. SUMMARY 

We summarize our detection of the surface brightness of 
EBL23 from resolved and unresolved galaxies fainter than 
V 555 = 23ABmag as follows (in ergs s ” 1 cm -2 sr ” 1 A -1 ): 

/ a (F 300W)= 4.0 (±2.5) xlO” 9 
/ a (F 555W)= 2.7 (±1.4) xlO” 9 
/ a (F 814W)= 2.2 (±1.0) xlO” 9 . 


The quoted errors are lcr combined uncertainties (statisti¬ 
cal and systematic) corresonding to 68 % confidence in¬ 
tervals, as described in §[y]. We can also define strict 
lower limits to the EBL from detected galaxies fainter 
than V 555 = 23 AB mag by computing the flux from all 
detected objects using the “aperture photometry” method 
described in §no|. The lower limits for EBL23 are (in ergs 
s ” 1 cm ” 2 sr”W4 -1 ): 

/23<v<27.5(F300W) > 3.2 (±0.22) xlO” 9 
/23 <v< 27 .5(F555W) > 0.89 (±0.01) xlO” 9 
l23<v<27.5(F814W) >0.76 (±0.01) xlO” 9 . 

For comparison with predictions of the EBL based on 
the local metal mass density and total star formation his¬ 
tory of the universe, the flux from galaxies brighter than 
V 555 = 23 AB mas should be added to these results, as 
shown in Figure |23|. We discuss these comparisons in Pa¬ 
per III. 
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APPENDIX 


MEAN ZODIACAL LIGHT FLUX THROUGH HST 
FILTERS 

The fiducial solar spectrum we have used for the pur¬ 
poses of creating a low-resolution spectrum of the zodiacal 
light from 2500A-l/im is a composite of the UV solar spec¬ 
trum of Woods etal. (1997), the optical spectrum from 
NL84, and the infrared spectrum produced by Arvesen 
etal. (1969), as recommended by Colina, Bohlin & Castelli 
(1996). The accuracy of the absolute flux of this solar 
spectrum is irrelevant to the accuracy of the zodiacal light 
spectrum, Izl(A); the absolute flux of /zl(A) is defined 
by the measured flux of the ZL in our own LCO spectra 
at 4650A, and by a combination of that measurement plus 
the color term, C(A, 4650A), at all other wavelengths. The 
color term itself is simply an empirical description of color 
of whatever fiducial solar spectrum we adopted relative to 
the observed color of the zodiacal light as measured in our 
FOS data. Thus, the absolute flux of /zl(A) is defined 
by the accuracy of the broad- band flux calibration of the 
FOS and LCO spectra, and the accuracy of the ZL mea¬ 
surement in the LCO spectra, as described in Paper II. 
The accuracy of the ZL measurement is, or course, depen¬ 
dent on the Solar Flux Atlas as discussed in Paper II. 

The spectrum Izl(A) as expressed in equation ^ is then 
an absolute flux-calibrated spectrum of the ZL, which we 
can convolve with the SYNPHOT throughput tables (us¬ 
ing the version released in May 1997) in the usual way 
to obtain the absolute flux of the ZL through each of the 
filters. The flux through the WFPC2 band is given by 


Iwf(A) 


/ T{\) I\(X) A dX 
f T(A) A dX 


(Al) 


in which T(A) is the effective throughput of a WFPC2 
filter (including telescope and detector efficiencies), I\ (A) 
is the spectrum of the zodiacal light, and all spectra are 
in units of ergs s ” 1 cm ” 2 sr ” 1 A” 1 . No additional error 
results from convolving the flux calibrated ZL with the 
bandpasses that define the WFPC2 system, as any error 
in the SYNPHOT synthetic photometry is incorporated 
in our estimate of the WFPC2 systematic uncertainty. 
See Paper II for a discussion of the LCO measurement 
of/ ZL (4650A). 


FLUX FROM THE LOW SURFACE BRIGHTNESS 
WINGS OF DETECTED GALAXIES 

In order to estimate the contribution from the wings of 
detected galaxies to the mean sky flux beyond 4ri so in a 
particular image, we have constructed a Monte Carlo sim¬ 
ulation. In this simulation, we sum the cumulative flux 
contributed by randomly placed, detectable galaxies to a 
given point on the sky, such that the total galaxy popula¬ 
tion simulated reproduces the appropriate surface number 
density as a function of magnitude. If a randomly placed 



26 


Detections of the Optical EBL: Results 



Fraction of mean sky flux 
0 0.005 0.01 



Fig. Bl.— Fitted values of scale length—^, and central surface 
brightnes s, u p, for galaxies from the Figures h9l (EBL data, open cir¬ 
cles) and pj (HDF data, points) assuming exponential light profiles. 
The HDF galaxies are plotted relative to the X-axis at the top of the 
plot; EBL galaxies are relative to the lower X-axis for clarity. The 
feature in both data sets which trails off to high h at constant fio 
corresponds to galaxies at the detection limit for which ^iso - /^O < 1 - 


galaxy is close enough to the “sky pixel” that the pixel 
would fall within the detection aperture of the galaxy, 
then the trial is rejected and another begins. The ex¬ 
tent of the galaxy apertures (how close galaxies can be 
to the sky pixel in question before that pixel falls within 
the galaxy’s detection aperture) are determined based on 
the mean “foreground” sky level and noise characteristics 
assumed for the image being simulated. The simulation 
continues until we obtain 10,000 “sky pixels.” 

As faint galaxies are weakly clustered on small projected 
scales (Colley etal. 1996, 1997; Roche etal. 1993), we have 
simply placed galaxies randomly in this simulation. For 
the surface density of galaxies as a function of apparent 
magnitude, we have adopted the HDF galaxy counts at 
V > 23AB mag, which go 1.5 mag fainter than our own 
data (see Figure |j). We have assumed that the observed 
light profiles of faint galaxies are adequately described by 
an exponential profile, p r = po exp (—r/h), where po is the 
central surface brightness and h is the scale length. Using 
the measured values of the core surface brightness in the 
central 4 (undrizzled) pixels, p CO re; the isophotal surface 
brightness, pj so , and the isophotal radius, ri so for the HDF 
galaxies, we can determine po and h (see Figures |l9] and 
20). As can be seen in Figures Bl and [TT^. h and po as a 


function of magnitude are consistent for our data and the 
HDF images, given the relative surface brightness limits 
of each. 


Fig. B2.— Histogram of flux par “sky pixel” found in two Monte 
Carlo simulations described in § |l0[ Each trial models the flux con¬ 
tributed to a random point on the sky by nearby galaxies. The lower 
X-axis shows flux in cgs units, while the top X-axis shows the flux 
as a fraction of the the mean foreground sky flux (dominated by ZL) 
at the levels observed in the EBL WFPC2 field. Surface brightness 
limits and sky noise levels used in this simulation to define the de¬ 
tection region around the simulated sources reflect the parameters of 
the V 555 EBL images. The thick line shows the sky pixel histogram 
which results from extending the “detection” aperture around sim¬ 
ulated source to the 1.4rj so . The thin line shows the histogram as¬ 
sociated with detection apertures extending to 4ri so . Vertical lines 
mark the mean “sky” flux identified in each simulation. 


gions associated with galaxies by the standard aperture 
size used for “total” magnitudes, y/2r -, so . In the second 
simulation, we extend the galaxy apertures to 4ri so . In 
the first simulation, we find that the flux from wings of 
galaxies contributes 2.3x 10” 10 ergs s -1 cm -2 sr -1 A -1 to 
the mean “sky” flux. As more galaxy light is excluded 
from the sky mean in the second simulation, the mean 
level of simulated “sky” pixels drops by just over 50% to 
l.lx 10 _lo ergs s -1 cm -2 sr -1 A -1 . The difference in flux 
in the “sky” pixels is within 10% of the flux we measured 
in the region 1.4 —4n so by ensemble aperture photometry, 
giving us confidence in our estimate of the flux coming 
from beyond 4ri so . Note that the total flux from beyond 
the standard galaxy photometry apertures, 2.3x 10 -10 ergs 
s^ 1 cm -2 sr -1 A -1 , is roughly < 0.5% of the mean sky 
level from diffuse zodiacal and galactic foregrounds, as in¬ 
dicated on the top x-axis of Figure [B2| . Note further that 
while the extragalactic pedestal from galaxies outside of 
the 4r; so detection apertures is only 0.1% of the ZL plus 
DGL foreground flux and 5% of the detected EBL in § pd[ 
it is 20% of the recovered flux from galaxies with V > 23 
AB mag by standard methods of galaxy photometry. 


EBL images 


In Figure B2, we plot histograms of the absolute flux 
contributed to 10,000 sky pixels at « 5500A in two simu¬ 
lations corresponding to the surface brightness limits of 
our F555W EBL images, m so = lcr s k y = 25.6V555 AB 
mag arcsec~ 2 . In the first simulation, we define the re¬ 


Comparison with the detection limits of the HDF 

As discussed in §[Tf], the crucial parameter for predict¬ 
ing the flux which will be recovered from a detected galaxy 
is the difference between the sky (limiting isophotal) sur¬ 
face brightness and the galaxy’s core surface brightness, 
A p = p, s ky — p-core- The surface brightness detection lim- 
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Fig. B3.— Each panel shows fitted h and po for galaxies in the 
HDF and EBL field s in unit magnitude bins. Axes and symbol types 
are as in Figure [Bl| . Lines show the /xo and h limits for galaxies with 
m = 23, 24, ...29. Galaxies in each magnitude bin can be seen to lie 
within the /xo, h relation for that magnitude range with the exception 
of galaxies near the surface brightness limits of the data, for which 
profile solutions become ill-defined. 


galaxies in the range 25 < I4o6 < 29.5 AB mag. This 
light will be included in the estimate of the local sky and 
thus subtracted from the flux within the galaxy aperture, 
doubling the total error. 

We have run the same Monte Carlo simulations de¬ 
scribed in §B1 using parameters that describe the noise 
and sky statistics of the HDF. Using fj,i so = l<j s k y = 
27.1U 6 o 6 AB mag arcsec -2 , we find an extragalactic contri¬ 
bution to the mean sky flux of 1 . 6 x 10 -10 and 0 . 8 x 10 _10 ergs 
s -1 cm -2 sr -1 A -1 for galaxy apertures extending to 
1.4ris 0 and 4rj so , respectively, as shown by the histograms 
in Figure B4. As for our EBL images, this is in good agree¬ 
ment with estimates of the fractional flux that should be 
recovered from galaxies as a function of magnitude for the 
corresponding values of A/j. In order to calculate the to¬ 
tal flux from galaxies detected in the HDF, we can sum 
the flux in the individual sources within standard, 1.4ri so 
apertures, and add to the resulting surface brightness twice 
the flux which lies outside of those apertures as identi¬ 
fied from the simulated sky histogram for HDF detection 
parameters. Doing so, we find that the total corrected 
flux from detected galaxies in the HDF 23 < 1406 < 31 
is 7.5x 10 ~ 10 ergs s -1 cm -2 sr -1 A -1 at 5850A. Convert¬ 
ing to the central wavelength of the V 555 filter, this is 
roughly 9.1xlO -10 ergs s -1 cm -2 sr -1 A -1 at 5250A, in 
good agreement with th e tot al flux from sources in our 


EBL field, estimated in §10.2 


Fraction of mean sky flux 
0 0.005 0.01 



Fig. B4.— Same as Figure B2 l but for surface bright detection 
limits, noise characteristics, and mean sky flux of the HDF V 606 
images. 


its of our data and HDF data are demonstrated in the 
Figures [if] and E§ Figure ||(] demonstrates that galaxies 
in the HDF with 1406 = 29 AB mag have A/r < l.Omag 
arcsec -2 on average, with some galaxies in that bin nearly 
reaching the limit p, s k y ~ /z C ore- If these galaxies have 
exponential profiles, then only ~ 20 % of their total flux 
can be recovered by standard photometry methods. By 
analogy with our data, galaxy apertures which extend to 
V2r iso in the HDF will miss at least 20% of the flux from 
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